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,  PART  T  W  0  -  T  HE  ,T  H  E  0  R  Y  AND 

PRACTICE  0F~  SOLVING  PR  OB  L  E  M  S  ,  I  N 
I  yN  T  E  R  N„A  L  BALLISTICS 

(T  H-/  0  RETICAL  AND  APPLIED 
P  Y  RODYNAK1  C  S) 

INTRODUCTION 

■  Oh  the  basis  of  the  widespread  study  of  the  phenomena  and 
processes  occurring  during  a ^discharge,  internal  ballistics  must. 


-  NS 


establish  the  laws  relating/the  conditions  of  loading  to  the 
quantifies  depending  upon  them  -  called  ballistic  elements  of 

discharge  — and  must  furnish  the  method  of  solving  a  large  number 

/ 

of  problems  encountered  in  practice. 

The  establishment  of  such  lavs,  providing  the  means  for  regulating 

/ 

a  discharge,  , constitutes  the  general  problem  of  internal  ballistics. 

The  conditions  of  loading  include  the  following:  the  dimensions 
of  the  powder  chamber  and  those  of  the  bore  of  the  barrel,  the  weight 
of  the  latter,  the  arrangement  of  the  rifling  in  the  bore,  the  weight 
rand  arrangement  of  the  projectile,  "the  pressure  necessary  to  overcome 
the  inertia  of  the  projectile,  the  weight  of  the  charge,  the  make  . 
of  powder,  the  physico-chemical  and  ballistic  characteristics  of 
the  powder,  the  characteristics  of  the  expansion  of  gases. 

The  ballistic  elements  of  a  discharge  include  the  path  of  the 
projectile  [  ,  its  velocity  v,  the  pressure  of  the  powder  gases  p, 
their  temperature  T,  all  values  varying  with  time,  and  also  the 
quantity  of  gasuf  formed  at  a  given  time.  .  . 

In  solving  the  above  general  problem  of  internal  ballistics, 
one  may  distinguish  two  fundamental  and  most  important  problems  of 
pyrodynaaics,  and  a  series  of  special  problems. 
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The  first  "fundamental  problem  consists  in  determining  by 
calculation  the  change  in  gas  pressure  and  the  velocity  of  the 
projectile  in  the  barrel  as  a  function  of  the  path  of  the  projectile 
and  of  time,  for  given  loading  conditions.  Together  with  the 
curves  p ,  / -v,  /- p,  t-y,  t  two  most  important  loading  characteristics 

1  of  the  gun  are  determined:  the  maximum  gas  pressure  p-  in-  the  bore, 

> 

and  the  muzzle  velocity  vA  of  the  projectile,  i.e.,  the  velocity 
of  the  projectile  at  the  instant  it  leaves  the  barrel  of  the  gun. 

This  problem  may  be  called  the  direct  problem  of  pyrodynaraics. 

For  given  conditions  of  loading  it  has  a  single  solution  -  a 
single  pressure  curve  with  maximum  pB,  a  fingle  velocity  curve  for 
the  projectile,  tad  a  muzzle  velocity  vA. 

By  varying  the  conditions  of  loading,  it  is  possible  to 
analyze  the  effect  of  these  conditions  on  the  variation  of  the  gas 
{  (  pressure  and  projectile  velocity  curves,  i.e.,  it  is  possible  to 

solve  a  series  of  special  problems  related  to  the  solution  of  the 
direct  problem. 

The  second  fundamental  problem  of  pyrodynamics  is  the  problem 
of  the  ballistic  design  of  the  gun;  it  consists  in  determining  the 
design  data  of  the  barrel  and  conditions  of  loading  necessary  to 
impart  some  definite  initial  (muzzle)  velocity  to  a  projectile  of 
a  given  caliber  and  weight.  This  velocity. is  determined  from  the 
tactical  and  technical  requirements  imposed  up<  n  the  gun  to  be 
t  constructed.  . 

* 

(  In  solving  such  a  problem,  the  maximum  gas  pressure  is  usually 

given. 

The  design  data  and  conditions  of  loading  insuring  that  a 
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projectile  of  a  given  caliber  and  weight  will  attain  the  desired 

Velocity,  are  obtained  from  the  solution  of  the  above  problem.  Once 

the  conditions  of  loading  are  given,  gas  pressure  and  projectile 

velocity  curves  are  drawn  as  a  function  of  path  and  time,  i.e., 

the  direct  problem  of  internal  ballistics  is  solved  for  the  selected 
■»  /  ,  • 
type  of  gun  and  charge. 

The  obtained  curve  p,  /  is  used  by-  the  engineers  to  calculate 
the  strength  of  the  gun  barrel  and  projectile  shell,  while  the 
curve  p,  t  is  used  to  design  the  carriage,,  the  time  fuzes  and  the 
Igniters.  At  the  same  time,  the  necessary  thickness  and  shape 
of,  the  powder  which  must  be  prepared  at  the  factory,  are  given. 

Thus  the  further  planning  of  the  entire  system  of  artillery 
and  of  the  necessary  ammunition  depends  to  a  considerable  extent 
upon  the  feasibility  and  rationality  of  the  selected  form  of  the 
ballistic  solution. 

This  is  why  the  problem  of  the  ballistic  design  of  guns  is 
the  principal  applied  problem  of  interior  ballistics; 

The  problem  of  ballistic  design  is  broader  than  the  first  problem; 
it  includes  the  latter  as  a  final  step  and  is  in  reality  an  inverse 
problem  of  interior  ballistics.  It  admits  of  numerous  solutions, 
numerous  combinations  of  gun  design  data  and  loading  conditions  under 
*  which  a  projectile  of  a  given  caliber  and  weight  will  attain  the 
required  muzzle  velocity. 

Because  of  the  indeterminate  character  of  the  solution,  there 
arises  the  need  of  developing  a  definite  method  for  obtaining  the 
necessary  answer  in  the  shorteot  possible  time,  and  for  selecting 
from  among  this  multiplicity  of  solutions  the  most  efficient  and 
desirable  solution,  satisfying  the  tactical  and  technical  requirements 
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imposed  upon  the  gun  to  be  designed. 

In  this  connection,  special  problems  arise  with  regard  to 
finding  the  most  desirable  solution,  and  for  obtaining  a  gun 

of  maximum  power  and  t-i nlmum  length  or  volume,  the  most  suitable 

,/ 

projectile,  and  the  most  desirable  load! ng,  conditions. 

The  solution  of  these  special  problems  permits  in  turn  to 
pose  the  problem  of  the  development  of  a  general  theory  and 
method  of  ballistic  design  which  would  take  into  account  the  most 
desirable  solutions  and  tactical  and  technical  requirements. 

Besides  the  indicated  fundamental  problems  of  internal 
ballistics,  there  is  also  a  series  of  special  and  secondary  problems 
introduced  below. 

For  a  given  bore  and  a  given  projectile  weight,  calculate  the 
weight  u>  of  the  projectile  insuring  a  given  muzzle  velocity  v^, 
and  the  thickness  2e^  of  powder  giving  the  required  maximum  pressure 

Pm*  .  _ 

Because  of  the  complexity  of  the  phenomenon  of-  discharge,  not 

all  of  its  details  can  be  taken  into  account,  even  approximately; 

some  of  these  details  must  be  neglected  and  can  not  be  introduced 

into  the  mathematical  equations  expressing  the  relations  between  the 

separate  processes  occurring  during  a:  discharge. 

For  this  reason,  the  equations  of  internal  ballistics  give  only 

* 

approximate  values  of  p>  v,  l ,  and  t.  But  since  in  practice 
these  equations  must  give  results  agreeing  with  experimental  data, 
it  is  necessary,  in  order  to  insure  this  agreement,  that  the  problem 
be  solved  by  selecting  certain  constant  characteristics.  When  these 
are  substituted  into  the  equations,  they  give  values  of  pB  and  vA 
F  -TS-7327-PJE  455 
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for  the  gases  and  the*  projectile,  respectively,  which  values  correspond 
to  the  results  of  firing  tests. 

The  yery  Banner  in  which  the  problem  is  posed  indicates  that 
the  processes  taking  place  during  a  discharge  are  not  yet  all 
known  arid  analyzed.  For  this  reason  one  of  the  main  problems  of 
\  internal  ballistics,  that  must  be  eventually  solved,  is  the  exact 

determination  of  constants,  those  of  the  gun  powder  in  particular, 
a®  derived  from  its  physical  and  chemical  properties.  The  determination 
of  the  powder  constants  involves  a  more  exact  method  of  pressure 
determination  by  experimental  means,  because  all  the  ballistic 
characteristics  (f,  a,  u^)  are  determined  from  the  latter. 

In  addition  to  the  problems  enumerated  above,  one  should  note 
the  problem  of  determining  the  variation  in  the  maximum  gas  pressure 
and  in  the  initial  projectile  velocity  under  specific  changes  in 
loading  conditions,  as  well  as  a  series  of  other  problems. 

The  fundamental  elements  of  a  discharge  -  /  ,  v,  p,  T,  and 

A  * 

t  -  are  interrelated  by  a  series  of  equations  expressing  the 
fundamental  processes  taking  place  during  a  discharge,-  i.e.  :  the 
burning  of  the  gun  powder  and  the  formation  of  gases,  the  trans¬ 
formation  of  the  thermal  energy  of  the  gases  into  the  kinetic 

\  *  - - — 

energy  of  the  system  projectile  -  charge  -  barrel,  and  the  movements 

of  parts  of  this  system. 

a 

The  methods  of  solution  of  theoretical  pyrodynamics  must  make 
.  it  possible  to  compute  and  establish  the  dependence  of  gas  pressure 
j  and  of  the  velocity  of  the  projectile  on  the  path  and  the  time  it 

takes  the  projectile  to  move  through  the  gun  barrel,  i.e.,  to 
solve  the  fundamental  direct  problem  of  internal  ballistics; 


F-TS-7327-RE 


456 


"A 


The  methods  of  solving  problems  in  pyrodynamics  may  bo  divided 
into  analytical/  numerical,  empirical  and  tabular  methods. 

i 

/The  empirical  methods  were  of. definite  advantage,  so  longas 
the  theoretical  concepts  of  internal  ballistics  had  not  been 
sufficiently  developed. 

They  were  bas'ed  on  some  relatively  simple  empirical  equations 
expressing  in  a  simplified  form  the  experimentally  obtained 

l 

interrelations  of  the  elements  of  a  discharge.  Tables  were  used 
along  with  these  equations,  which  tables  offered  the  means  for 
computing  very  rapidly  the  elements  of  the  curves  depicting  gas 
pressure  and  projectile  velocity. 

The  empirical  methods  were  derived  from  the  analysis  of 
experimental  data  obtained  in  firing  weapons  under  different 
conditions,  with  the  characteristics  and  constants  entering  into 
these  expressions  determined  from  the  conditions  of  the  experiment. 

The  disadvantage  of  these  methods  (formulas  and  tables)  consists 
in  the  fact  that  they  fail  to  take  into  account  certain  very  important' 
factors  and  conditions  of  loading,  and  that  such  methods  may  be 
applied  only  under  the  conditions  and  within  the  limits  established 
for  the  given  case. 

The  number  of  empirical  equations  and  tables  is  very  large; 
prior  to  the  development  of  analytical  solutions,  they  were  of 
primary  value  because  of  their  simplicity.  But  the  appearance  of 
exact  theoretical  solutions,  taking  into  account  with  sufficient 
completeness  the  influence  of  most  of  the  conditions  of  loading  and 
singularities  of  the  processes  occurring  during  a  discharge,  made 
it  possible  to  solve  all  the  fundamental  problems  of  pyrodynamics 
by  means  of  exact  analytical  relations.  As  a  result,  many  empirical 
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‘equations  and  tables  have  lost  their  significance  and  are  now  used 
only  in  certain  auxiliary  cases. 

The  analytical  methods  are  based  on  a  series  of  assumptions 
characterizing  the  conditions  of  powder  burning  and  the  motions 
of  the  gases,  projectile  and  gun;  these  assumptions  are  based  chiefly 
on  experimental  or  theoretical  data  expressing  the  physical  side 
of  the  process  of  discharge. 

For  this  reason,  analytical  methods  of  solution  give  a  more 
profound  understanding  of  the  real  n£'  ure  of  the  phenomenon  than 
empirical  methods,  and  approach  more  closely  the  essence  of  the 
processes  taking  place  during  a  discharge. 

In  the  analytical  method  the  problem  is  reduced  to  the-  solution 
and  integration  of  differential  equations  of  different  types.  This 
solution  can  be  obtained  with  greater  accuracy  (in  which  case  the 
resulting  equations  become  more  complex),  or  approximately  (which 
results  in  simpler  relations). 

I 

Solutions  may  be  given  for  the  more  complex  cases  obtained  in 
practice,  and  also  for  simplified,  admittedly  schematic  cases,  in  which 
case  the  analysis  of  the  relationships  is  simplified. 

Solutions  may  be  based  on  the  geometric  and  physical  laws  of 

\ 

powder  burning.  .  ^ 

Tables  of  auxiliary  values  or  functions,  necessary  to  calculate  1 
certain  intermediate  values,  are  prepared  in  order  to  expedite  and  \ 

simplify  the  computations  involved  in  the  solution  of  problems.  \ 

\ 

Numerical  methods  of  integrating  a  system  of  differential 
equations  are  used  along  with  the  analytical  methods^.  The  integration 
is  usually  performed  by  the  method  of  finite  differences  or  by 
expansion  in  Taylor's  series.  These  methods  are  resorted  to  in 
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especially  difficult  cases,  when  the  value  of  one  or  several  parameters 
varies  throughout  the  process  of  discharge  and  their  variation  does 
not  permit  to  solve  the  problem  analytically  in  finite  form.  This 
happens,  for  example,  when  the  cross  section  of  a  barrel  bore  varies 
.  (tapered  bore),  or  when  the  parameter  6  varies  throughout  the 
discharge  accompanied  by  a  varying  gas  temperature  or  by  a  change 
of  the  coefficient  depicting  to  secondary  work  done  in  the  process, 
etc. 

On  the  basis  of  analytical  or  numerical  solutions,  it  is  possible 
to  set  up  numerical  tables  of  the  fundamental  elements  (p,  v,  l  ,  t) 
for  different  loading  conditions  and  some  general  constants.  These 
tables  enable  one  to  plot  very  rapidly  the  necessary  curves  p,  /  and 
v,  !  or  p,  t  and  v,  t  with  a  minimum  number  of  calculations.  In  so 
doing,  the  process  of  solving  the  direct  problem  is  greatly  simplified 
and  expedited.  These  tables  are  usually  set  up  for  certain  average 
values  of  the  constants  (characteristics  and  shape  of  the  powder), 
although  in  practice  one  may  encounter  a  series  of  regressions  from 
these  average  values.  In  that  case  it  is  necessary  to  introduce 
appropriate  corrections  into  the  results  obtained. 

Thus  the  ballistic  tables  for  the  solution  of  direct  problems 
of  pyrodynamics  are  in  reality  analytical  equations  reduced  to 
numerical  values  in  a  series  of  concrete  loading  conditions. 

However,  the  ballistic  tables  enable  one  to  solve  a  series  of 
problems  which  cannot  be  solved  directly  by  means  of  analytical 
equations. 

The  basic  difference  between  tabular  values  and  analytical 
equations  is  the  following:  because  of  their  complexity,  the 
analytical  equations  do  not  give  a  direct  relationship  between 
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pressure  or  velocity  and  path  length,  for  example;  these,  variables 
are  usually  related  through  some  auxiliary  variable. 

In  tables,  on  the  contrary,  the  basic  elements  of  discharge 
are  interrelated  directly:  the  pressure,  the  projectile  velocity, 
and  the  time  of  its  travel  through  the  barrel  are  given  .in  function 
of  the  path  traversed  by  the  projectile;  this  simplifies  considerably 
the  analysis  and  permits  the  development  of  a  special  method  for 
solving  problems  which  cannot  be  solved  by  analytical  means. 

The  development  of  the  theory  of  ballistic  design  became 
possible  only  with  the  introduction  of  tables  for  the  solution  of 
internal  ballistic  problems. 

For  this  reason  the  first  tables  prepared  in  our  country  by 
Prof.  N.F.  Drozdov  on  the  basis  of  his-  exact  solution  given  in  1910, 
are  of  great  importance.  These  very  tables  simplified  and  expedited 
the  calculations  involved  in  the  ballistic  design  of  weapons,  and 
gave  the  engineers  a  reliable  means  of  solving  rapidly  inverse 
problems  in  pyrodynamics. 

They  also  served  as  an  example  for  a  series  of  more  detailed 
tables  compiled  subsequently. 


( 
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SECTION  SIX  -  A  N  A  L  Y  TXI  C  A  L 
METHODS  OF  SOLUTION  OF  THE 
DIRECT  PROBLEM  OF  INTERNAL 

BALLISTICS. 

BASIC  ASSUMPTIONS 

When  we  examined  the  phenomenon  of  a  discharge,  we  had  pointed 
out  its  extreme  complexity  and  the  fact  that  some  of  the  factors 
influencing  the  results  were  still  insufficiently  known.  For  this 
reason,  when  solving  theoretically  the  fundamental  equation  and 
deriving  the  relations  between  the  physico-chemical  and  mechanical 
phenomena  in  a  discharge,  it  is  necessary  to  take  recourse  to  certain 
simplifications  and  schemes. 

The  basic  assumptions  are  as  follows: 

1)  The  burning  of  powder  obeys  the  geometrical  law  of  combustion. 

2)  The  powder  burns  under  an  average  pressure  p. 

3)  The  composition  of  the  products  of  combustion  does  not  change 
during  burning,  nor  during,  the  adiabatic  expansion  of  the  gases  (f  - 
and  a  are  constant)  after  the  powder  is  burned, 

4)  The  rate  of  burning  is  proportional  to  the  pressure: 


u  -  u1p. 

5)  The  auxiliary  work  done  is  proportional  to  the  principal 
work  of  '.he  forward  motion  of  the  projectile,  and  is  represented 
by  the  coefficient^. 

6)  The  projectile  starts  moving  when  the  pressure  developed  in 
the  chamber  by  the  partial  burning  of  the  charge  equals  p^,  i.e., 
when  the  pressure  is  sufficient  to  force  -the  driving  hand  completely 
into  the  rifling  of  the  bore;  the  gradual  forcing  of  the  band  and  the 

t 

increasing  resistance  encountered  by  it  are  not  taken  into  account. 


ss 
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7)  The  work  done  in  forcing  the  driving  band  is  not  accounted 

* 

for  separately,  nor  the  increasing  velocity  of  the  projectile  during 
the  gradual  forcing  of  the  band. 

8)  The  expansion  of  the  barrel  during  the  discharge,  the 
gases  escaping  through  the  clearance  between  the  driving  band  and 
the  wajls  of  the  gun,  and  the  air  resistance  are  disregarded. 

9)  The  cooling  of  the  gases  through  heat  transfer  to  the  walls 
of  tne  barrel  is  no.t  accounted  for  directly,  and  may  be  taken  into 
account  indirectly  (for  example,  by  decreasing  the  force  f  -  RT^ 

or  increasing  ©  -  l/(A  +  BTav<). 

10)  The  motion  of  the  projectile  is  considered  only  until 
it  passes  the  muzzle  face. 

11)  The  quantity  ©  -  (c^/c^)  -  1  is  taken  as  its  average  value, 
constant  throughout  the  discharge. 

The  assumptions  enumerated  above  make  our  representation  of 
a  discharge  more  schematic,  and  deviate  the  phenomenon  to  a  greater 
or  lesser  degree  from  reality.  For  this  reason  the  relations  obtained 
in  the  solution  will  express  the  physico-mechanical  nature  of  the 
discharge  only  with  a  certain  degree  of  approximation.  Thus  the 

values  of  the  fundamental  elements  (maximum  gas  pressure  p  ,  initial 

m 

or  muzzle  velocity  vA  of  the  projectile)  obtained  from  these  equations 
may  not  coincide  with  the  values  obtained  experimentally.  Nevertheless, 
in  order  to  solve  practical  problems,  it  is  necessary  to  obtain 
analytical  data  which  would  agree  with  experiment.  For  this  reason 
(keeping  in -mind  the  complexity  of  the  discharge  phenomenon,  the 
incomplete  knowledge  of  its  elements,  and  the  disagreement  between 
our  basic  assumptions  and  reality)  it  is  necessary  to  introduce 
coefficients  of  "agreement  with  experiment"  into  the  constants  obtained. 
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Such  a  method  is  widely  used  in  various  scientific  laws 
(hydrodynamics,  aerodynamics,  etc.)  dealing  with  complex  phenomena, 
whose  details  cannot  be  fully  analyzed. 

Eventually,  as  our  knowledge  is  further  developed,  we  will  find 
it  possible  to  render  some  of  the  assumptions  with  greater  accuracy 
and  take  into  account  some  of  the  conditions  not  yet  understood. 

As  new  experimental  data  is  accumulated  arid  new  methods  are  applied, 
the  deductions  arrived  at  may  be  modified  and  even  replaced  by 
others  of  a  more  complete  and  exact  nature. 

In  solving  the  fundamental  equation  of  pyrodynamics,  one  should 
strive  to  obtain  the  maximum  possible  mathematical  accuracy.  However, 
in  that  case  some  of  the  expressions  become  excessively  cumbersome, 
so  that  even  exact  formulas  will  fail  to  represent  the  true  phenomena 
of  a  discharge;  and  for  this  reason  certain  simplifications  may  be 
used  with  advantage  in  the  process  of  solution. 

A  comparison  of  these  simplified  solutions  with  the  exact  ones 
may  show  the  extent  of  mathematical  error  involved  with  the  use 
of  the  same  constants  and  conditions. 

With  an  appropriate  selection  of  constants,  the  somewhat 
simplified  solutions  may  also  yield  results  approaching  experimental 
data  as  closely  as  those  obtained  by  the  use  of  more  exact  equations. 

CHAPTER  1  -  SOLUTION  OF  THE  FUNDAMENTAL  PROBLEM  WHEN  THE 
PRESSURE  TO  OVERCOME  THE  PROJECTILE  INERTIA  IS  KNOWN,  AND  WHEN 
BURNING  PROCEEDS  ACCORDING  TO  THE  GEOMETRIC  LAW. 

As  we  have  shown  above,  the  fundamental  equation  of  pyrodynamics 
includes  a  large  number  of  constants  characterizing  the  projectile, 
charge  and  powder  which  determine  the  conditions  of  loading,  and  the 
four  variables,  y,  v,  l  and  p,  which  are  called  the  elements  of  a  shot 
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In  order  to  establish  the  relation  between  the  elements  of  a 
shot,  new  equations  are  added  to  the  fundamental  equation,  which 
•re  the  equations  of  powder  burning  and  projectile  motion;  this 
leads  to  the  appearance  of  a  new  variable,  the  time  t,  and  to  the 
appearance  of  quantity  z  when  burning  proceeds  according  to  the 
geometrical  law. 

We  obtain  as  a  result  the  following  system  of  equations: 

The  fundamental  equation  of  pyrodynamies : 

8 


ps(^  +  l)  -  fu>*f  -  — <f mv  . 

2 


(1) 


The  rate  of  powder  burning: 


u 


de 

dt 


-  uxp. 


(2) 


*  ( 


The  law  of  generation  (inflow)  of  gases: 


«  Y.  z(l  +  Xzj  «  x  z  +  xXz  . 


(3) 


The  law  of  motion  of  the  projectile: 


ps  -  <fm 


dv 

dt 


(4) 


or 


ps  *  if  mv 


dv 

dt 


(5) 


The  totality  of  these  equations  affords  the  solution  of  the 
fundamental  mathematical  problem:  of  determining  the  curves  p, [ 
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and  v,  l  and  also  p,  t  and  v,  t  and  of  finding  in  particular  the 
maximum  gas  pressure  pm  and  the  muzzle  velocity  of  the  projectile. 

We  first  solve  the  problem  for  regressive  powder  shapes, 
using  the  two-term  formulas  (X  >  1,  X  <  0,  p.  =  0)  and  the  assumptions 
enumerated  above. 

We  shall  solve  the  problem  for  all  the  periods  of  a  shot  in 
succession. 

1.  PRELIMINARY  PERIOD 

When  establishing  the  relations  for  this  period,  we  shall  assume 
the  simplest  form  of  the  phenomenon:  the  instantaneous  forcing  of 
the  projectile  band  into  the  rifling. 

Fundamental  assumption.  If  the  force  necessary  to  overcome  the 
resistance  encountered .by  the  driving  band  of  the  projectile  in 
completely  penetrating  the  rifling  is  and  the  cross  section 

of  the  bore  is  s,  the  quantity  pQ  -  IIq/s  will  be  called  "the  pressure 
to  overcome  the  inertia  of  the  projectile"  or  the  forcing  pressure. 

We  shall  assume  that  the  projectile  is  set  in  motion  at  tl\e  instant 
the  gas  pressure  attains  the  value  p^. 

Up  to  that  moment  the  burning  of  the  powder  takes  place  in 
a  constant  volume.  For  this  reason  the  preliminary  period  may  be 
called  pyrostatic  and  one  may  apply  the  already  known  equations  of 
pyrostatics. 

In  this  period,  besides  the  forcing  pressure  p^,  we  will  be 
interested  in  the  portion  of  the  charge  burned  at  the  instant  the 
projectile  i^  set  in  motion,  in  the  relative  thickness  of  the  powder 
Zq  -  eQ/e^,  and  in  the  relative  surface  area  of  the  powder  Sq/S^  -  6^. 

These  quantities,  characterizing  the  end  of  the  preliminary 
period,  are  simultaneously  the  initial  values  of  the  first  period. 
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Let  us  introduce  the  fundamental  equation^  for  the  preliminary 


period. 


The  igniter  is  burned  first  and  the  pressure  developed  in  the 

/ 

chamber  is  Pg,  which  may  be  computed  by  the  following  formula: 


PB  " 


fBwB 


If  ~  _  >a„un 

"0  b  B  B 


where  Wq  is  the  volume  of  the  chamber;  w/3  is  the  volume  of  the 
charge  proper;  and  fg,  afi, w ^  are,  respectively,  the  force,  co¬ 
volume,  and  weight  of  the  igniter.  Under  the  usual  conditions  of 
ignition,  cigUg  may  be  neglected. 

The  charge  proper  will  ignite  when  the  pressure  reaches  p  ; 

O 

at  this  instant  the  pressure  is  determined  by  the  general  equation 
of  pyrostatics,  which  takes  into  account  the  effect  of  the  igniter. 
At  the  instant  the  driving  band  is  forced  in  the  rifling,  a  certain 
portion  of  the  charge  will  have  burned,  and: 


po  -  ■’b  + 


V- 


u  u 


h  6  j^O 


where 


1  _  1 
Tq  "  °  “  t- 


Inasmuch  as  the  forcing  pressure  p^  i«,known(^),  we  can  determine 

(*)  Po  varies  between  250  and  400  kg/c»2  for  shells  and  between 
300-500  kg/cm2  for  bullets  when  the  entire  side  surface  is  forced  into 
the  rifling  of  the  bore. 
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what  part  ^ -of  the  charge  will  have  been  burned  at  the  instant  the 
projectile  is  set  in  motion.  Solving  equation  (7)  for  v^»  ,  we  obtain: 


/ 

/ 

If  we  may  neglect  the  pressure  of  the  igniter,  because  pQ  is 
known  only  approximately,  while  is  small;  we  will  obtain  a 
simpler  expression  for  computation: 


-The  quantity  v^q  mainly  depends  on  A  and  varies,  in  general, 
between  0.02  and  0.10. 

If  the  amount  of  the  burned  portion  of  the  powder  is  known, 

and  the  law  of  powder  burning  is  in  the  form: 


'I'  "X  z(l  +  Xz)  -X  z  +  xXz^ 


we  can  determine  the  relative  thickness  z -  e^/e.^  of  the  powder 
burned  at  the  start  of  motion,  and  the  relative  surface  area  6q. 
We  find  6 q  from  the  following  formula: 
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r.U  ;•>:  -  '* 


««*  3J 1  C'i-V_  >  ^  P*  „  ->,  J  -jC* 


and  zn,  from  the  equation  G>n  -  1  +  2Xzft: 


s  -  1  (60  -  1)(S  +i)  ej  -  1  2+ 

X  «  - - - —  “  - - - - -  “  -  -  - - - , 

U  2X  2X((50  +  1)  2X(6q  +  1)  (G0  +  l)x 


Since  usually  S^sssl,  the  following  approximation  is  correct 


Besides  these  characteristics,  we  will  also  require  the  value 
-  the  length  of  the  free  space  in  the  chamber  at  the  start  of 
motion.  This  reduced  value  is  determined  from  one  of  the  following 


expressions : 


Lo  „  ± 


G  "0 


“  _'2° ).  1  ±  .  Jl* 

l  hi  0  \  *  *i  0 


where 


-  I 


4  U  -  ±  -  is  V 

0  \a  .  S  h }’ 


I  U  W  co 

~ -  -  -  — . 

0  8  Wq  s 


2.  FIRST  PERIOD 


In  deriving  the  fundamental  relationships  for  ithe  first  period. 

Prof.  N.F.  Drozdov  was  the  first— to  propoee  the  introduction  of 

\ 

a  new  independent  variable,  x  -  z  -  z q  (the  relative  thickness  of 
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the  powder  burned  after  the  projectile  is  set  in  motion). 

At  the  instant  the  projectile  is  set  in  motion  z  -  zQ  and 

x  «  0;  at  the  ending  of  burning  z  -  1  and  x  »  1  -  zn. 

K  K  d 

Thus  the  limits  of  variation  of  the  new  argument  are  known  in 
advance.  Let  us  express  all  four  fundamental  elements,^.,  v,  /  , 
and  p,  as  a  function  of  this  argument. 

Relation  'f' “  f^(x).  Substituting  z  -  z^  +  x  in  the  formula 
S'  “  *  z  ••  xXz^,  we  obtain: 

S'  "  K  z0  +  K^zo  +  +  2Xz0)x  +  X  Ax2, 

but 

*z0  +  xAz*  .  to;  !  +  2Xz0  -  G0. 

Introducing,  according  to  Drozdov,  the  additional  designation 
*G0  "  kl*  we  obtain  the  desired  relation: 

S'  -  +  k1x  +  kAx2.  (10) 


2.  Relation  v  -  f2(x).  The  velocity  v  enters  into  the  equation 
of  motion: 


dv 

Sp  -  __ 


In  order  to  eliminate  p  and  t,  we  add  the  law  governing  the 


rate  of  burning: 
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Multiplying  these  equations  term  by  term  and  simplifying, 
we  obtain: 


dv 


s  de 

y<?a  el 


sc, 

x 

<f  5iU  ^ 


dz 


si, 


m 


dz. 


Integrating  from  0  to  v  and  from  Zq  to  z: 


<f  M 


— x. 
«?  m 


Prior  to  the  end  of  burning 


(11) 


v  « 
K 


si, 

cp  m 


•d  -  V 


K 

J.  _5_  (1  -  z„). 
q  :  s  0 


(12) 


Consequently,  the  velocity  of  the  projectile  at  the  end  of 
burning  can  be  computed  in  advance,  if  the  impulse  of  the  powder 
pressure  I %  *■*  e^/uj  and  the  cross-sectional  loading  of  the 
projectile  are  known. 

Inasmuch  as  the  quantities  and  1  -  Zq  vary  relatively  little, 
the  velocity  of  the  projectile  at  the  end  of  burning  of  the  powder 
depends  in  the  main  on  the  ratio  of  the  impulse  I  to  the  cross- 

A 

sectional  load  q/s  on  the  projectile,  and  during  burning,  the  velocity 
of  the  projectile  varies  in  proportion  to  x. 

Equation  (12)  permits  to  compute  v^t  but  it  does  not  tell  us 
the  point  on  the  path  the  projectile  at  which  the  powder  is  burned, 
whether  the  speed  v^  is  properly  chosen  for  the  given  gun,  or 
whether  the  powder  is  fully  burned  before  the  projectile  leaves  the 
gun.  For  this  reason,  this  equation  alone  is  insufficient,  and  it  is 
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necessary  to  find  also  the  equation  for  the  path  traversed  by  the 

*  . 

projectile  at  the  end  of  burning. 

Equation  (12)  is  plotted  in  fig.  133;  it  shows  the  curve 
and  gives  the  value  of  v  ,  but  it  does  not  show  the  position  of 
the  projectile  at  the  end  of  burning. 


Fig.  133  -  Path  of  the  Projectile  at  the  End  of  Powder  Burning. 

3.  Relation  l  -  /^(x).  In  order  to  determine  the  path  of  the 
projectile,  two  equations  must  be  used:  the  fundamental  equation 
of  pyrodynamics  and  the  equation  of  motion  of  the  projectile  in 
the  form  of  elementary  work: 


psd/  **cfmvdv. 


2 

where  V|j  *»  2fu/p0m.(*) 

We  eliminate  p  by  dividing  the  second  equation  by  the  first; 


(*)  The  Russian  subscripts  D p  denote:  path  traversed  inside  a  barrel. 
Editor.  f 
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Since  v  and  ^  are  functions  of  x  according  to  equations  (10) 
and  (11),  the  right-hand  side  of  this  differential  equation  may  be 
represented  as  a  function  of  x.  . 

Designating  this  function  by  dF(x)  and  substituting  for  v  and  ^ 
their  expressions  in  x  we  obtain: 


dF(x) 


<f  2m2 


2  2 

2  s  fg'fme 

t0  +  k  x  +xXx  -  — - — -x 


s2I2 
8  aK 

f  wcf  m 


t0  +  V  " 


S2I2 
1K  8^ 

fujtfrn  2 


2  o  / 

The  same  group  s  I^/fwym  of  constants  and  characteristics 
appears  in  the  numerator  and  the  denominator.  As  suggested  by  Prof. 

N.F.  Drozdov,  it  is  represented  by  B  and  is  called  "the  parameter 

/ 

of  the  loading  conditions"  (Prof.  N.F.  Drozdov's  parameter): 
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where : 


i  (x)  -  x2  -  :~x 
1  B1 


Prof.  N.F.  Drozdov  was  the  first  to  solve, this  equation  exactly, 
in  1903,  by  reducing  it  to  the  form  of  a  linear  equation  of  the 
first  order: 


d  (  ^  B  *  |  ^  _ B  x 

dx  Bj.  Ei(x)  "  '  B  Ijtn) 


+  V  -Qx, 


where  P  and  Q  are  functions  of  x. 

X  X 


The  full  solution  of  this  equation  is  presented  later. 

A  simpler,  but  approximate  solution  is  obtained  if  we  assume 
If  -  "  const. 

It  will  be  presented  later  with  the  designation  of  the  parameters, 
with  some  of  the  auxiliary  functions  derived  according  to  Prof. 

Drozdov. 

During  burning  of  the  powder  at  the  start  of  the  projectile’s 

motion,  l ^  varies  within  the  limits  of  t  and 
T  To 

lu>li>h’ 
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where 


The  rate  of-^'change  of 
Assuming  that 


increases  with 

“  l  and  integrating  equation  (13’)  we  have: 
•  av. 


i 

f 

j- 

0 


dl 


xdx 


The  integral  of  the  right-hand  side  is  obtained  by  decomposing 
the  integrand  into  the  simplest  fractions;  it  is  a  logarithmic 
function  of  x  which  we  shall  temporarily  designate  by  In  Z^.  The 
left-hand  side  is  integrated  also: 

l  \  B 

1  +  — : — - In  Z  , 

/  /  Bl 

Yav/ 

whence : 

B 

l  (zx  1  -1).  (14) 

Thus  the  expression  for  the  path  l  as  a  function  of  x  is  more 
complex  than  the  expressions  for  and  v. 

Substituting  into  it  x  -  1  -  z^,  we  can  find  the  path  at 
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the  end  of  powder  burning.  Comparing  it  with  the  full  path  lA 

s 

traversed  by  the  projectile  within  the  bore,  it  is  possible  to 
determine  whether  the  thickness  of  the  powder  and  the  velocity 
of  the  projectile  are  correctly  chosen  for  the  given  gun. 


In  computing  l...  we  may  use  in  the  equation  t 


av. 


av. 


"  l 


* 

\  1 

l:  -  A 

fa  -A 

J  'av. 

s 

\  1 

the  expression 


av. 


-  yo  +T 


/ 


'*•  »»»»'«?  "'te  r+*. 


The  investigations  of  Prof.  G.V.  Oppokov  in  his  book  ”0 

TOCHNOSTI  NEKOTORYKH  ANALITICHESKIKH  SPOSOBOY  RESHENIYA  OSNOVNOI 

x  / 

ZADACHI  VNUTRENNEI  BALLISTIKI  DLIA  PERVOGO  PERIODA"  (Concerning 

the  Accuracy  of  Certain  Analytical  Methods  of  Solving  the 

Fundamental  Problem  of  Internal ‘Ballistics  for  the  First  Period), 

2  7 

1932^  ~  have  shown  the  following.  When  the  loading  density  is 

A  -  0.5-0. 7,  formula  (14)  is  very  accurate  for  evaluating  pB  and  v^, 

if  l.  is  not  taken  to  have  the  same  value  for  all  the  values  x 
'  av. 

from  0  to  1  -  Zq,  and  if  a  different  value  of  /  is  taken  for 

'  av. 

every  value  of  x,  assuming  either  of  the  following  values  for  L, 

'av. 

in  the  formula:  -  (t  -  H'q)/2  (Oppokov)  or  ^av<  -  + 'f*)/2 

(Serebryakov) c 

Inasmuch  as~r  ls  directly  proportional  to  v  /"equation  (11  )_7> 
equation  (14)  gives  in  fact  the  direct  relation  between  the  path  / 
and  the  projectile  velocity  v. 

The  expression  for  Z%  presented  below  shows  that  this  relation 
is  expressed  by  a  rather  complex  function. 

4,  delation  p  -  f^(x).  The  pressure  p  is  found  from  the 

fundamental  equation  of  pyrodynamics :  ■ 

1 
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p 


fw 

S 


V  1 


If  the  quantities^,  v,  and  l  are  replaced  in  the  right-hand  Side 
'by  their  expressions  ir.  function  of  x,  then 


P  - 


fu>  to  *  ~  Bxx 

~s  "  "  *  B 

'  +  L.  (z*  Bi 


f 

s 


'"X 

1  av. 


-  1) 


.  B6  2 

d,  _  - x 

2 


+  / 


(15) 


where  ty  can  be  represented  as  a  function  of  x  as  well. 

Inasmuch  as^,  v,  and  l  are  already  determined,  it  is  no 
longer  necessary  in  computing  p  to  use  equation  (15)  which  is 
expressed  in  terms  of  x,  and  the  numerical  values  of  v|>,  v,  \  , 

and  l  can  be  substituted  in  the  preceding  equation. 

By  attributing  to  x  different  values  within  the  limits  of  0 
to  i  -  Zq,  equations  (10),  (11),  (14),  and  (15)  permit  one  to 
f  ind  the  values  of  all  the  elements  y,  l  ,  and  p  of  a  shot  entering 
into  the  fundamental  equation  of  pyrodynamics,  and  to  plot  a  curve 
showing  the  variation  of  p,  v,  4^  as  a  function  of  /  ,  i.e.,  the  curves 
p,  /  and  v,  /  . 

Consequently,  the  proposed  problem  concerning  the  solution  of 

the  fundamental  equation  of  pyrodynamics  has  been  resolved,  and  the 

relation  between  the  elements  has  been  found. 

*  * 

Substituting  the  value  of  xR  «  1  -  In  the  above  equations,  we 
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find  all  the  elements  v^.,  p^f  and  l  corresponding  to  the  instant  the 
burning  of  the  powder  ends  (4* -  I)*  these  values  will  be  the  initial 
values  in  the  second  period. 

Note.  The  expression  for  the  projectile  velocity  may  be 
replaced  by  the  following: 


2t2  A2 

2  E  lk  o  s  *K  b  2fw  2  BO  2  2 

ii»,»  ‘f"'8  2  np 


whence. 


This  expression  brings  out  the  effect  of  the  limiting  velocity 

of  the  projectile  and  that  of  the  parameter, of  loading  conditions,  B. 

Since  in  most  guns  B  varies  within  narrow  limits,  it  follows  that 

the  velocity  mainly  depends  upon  the  potential  f/e  of  the  powder 

and  upon  the  relative  weight  >/q  of  the  charge. 

x 


y  xdx/^(x) 


Determination  of  the  function  Z  re1 
— - - * 

evaluate  the  integral  \ 


~Xd2L_>a 

5i(x) 


In  order  to 
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but: 


and/  consequently: 


b  +  1 


and  x2  are  expressed' by  equations  (16). 

Substituting  .here  these  values  of  x1  and  Xj>,  we  get: 

.  b+1  b-1 


Inasmuch  as  the  quantity  b  -  y  1  +  4ByVy/k^  “  y  1  +  4^ 

is  itself  a  function  of  the  parameter  y  -  the  function 

Zx  actually  depends  only  upon  two  quantities:  the  constant 
2 

-  B^g/k^  and  the  variable  0  -  B^X/k^. 

From  these  data  it  is  possible  to  set  up  a  table.  Since  the 

-B/Bi 

equation  of  the  path  contains  the  expression  Z„  the  tables  ai*e 

X 

set  up  for  log  Z  *  to  make  their  use  more  convenient. 

The  quantities  entered  (introduced)  are  y  and  p. 


It  is  not  difficult  to  show  by  another  method  that 


l  2  *1  ^0 

Ox-  — — x—  — 
B1 


A 

-  In  Zx  is  a  function  of  y  -  and  P  **  B^x/k^,  if  the  numerator 

2/9 

and  denominator  of  the  integrand  are  multiplied  by  Bj/k^.  Then, 


This  expression  actually  shows  that  In  Zx  is  a  function  of  y 

and 

-1 

The  table  of  the  logarithms  of  the  function  (log  Z^  )  is 
presented  below  (Table  1). 
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Procedure  for  Using  the  Table. 

For  every  problem  we  will  have  one  value  for  the  entry  parameter 
2 

f  "  BlVkl  and  a  series  of  values  p  -  BjX/kj,  where  x  varies 

between  0  and  1  -  %  . 

0 

When  determining  log  Z~\  write  down  the  values  from  the  -columns 
containing  the  nearest  smaller  and  larger  tabular  values  of  y, 
so-  that  the  value  of  y  obtained  from  the  solution  would  fall  between 
them,  the  coefficient  of  interpolation  will  be  the  same  along  all 
the  horizontal  rows.  For  this  reason,  it  is  more  convenient  to 
interpolate  first  along  the  horizontal  rows  between  which  are  contained 
the  values  of  p  selected  in  the  problem,  and  then  to  interpolate  along 
the  columns  (vertically)  using  the  corresponding  interpolation 
coefficients  P. 

In  order  to  reduce  the  number  of  vertical  interpolations  (except 

such  cases  when  log  Z”1  is  used  for  computing  the  values  p  and  p  ), 

®  K. 

it  is  more  convenient  to  assign  tabular  values  of  p  for  the 

intermediate  values  of  x  and  to  perform  only  the  horizontal  interpolation 

for  y.,  and  then  determine  x  by  means  of  equation  x  «  k^P/B 

3.  DETERMINATION  OF  THE  MAXIMUM  PRESSURE  GENERATED 
BY  THE  POWDER  GASES 

The  maximum  gas  pressure  pm  in  the  barrel  is  the  most  important 

ballistic  characteristic  of  a  gun.  Its  value  depends  on  the  chosen 

conditions  of  loading,  and  the  obtainment  of  the  desired  value  of  d 

^m 

serves  as  a  criterion  or  control  for  the  proper  selection  of  the 

weight  of  the  charge,  the  thickness  of  the  powder,  and  other  loading 
conditions. 

For  this  reason,  it  is  sometimes  important  to  be  able  to  compute 
the  pressure  pra  for  the  given  loading  conditions,  without  constructing 
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the  entire  p,  (  pressure  curve.  In  order  to  achieve  this,  it  is 
necessary  to  derive  first  a  formula  for  determining  the  value  of 
x  for  which  the  gas  pressure  is  maximum. 

In  this  case  the  derivative  dp/d(  or  dp/dt  must  be  equated 
to  zero.  The  expression  for  the  derivative  was  derived  earlier 
by  differentiating  the  expression  for  p  from  the  fundamental 
equation  of  pyrodynamics. 


dP 

dl 


P  f iu  X  6 
+  /  t  s  lK  v 


1  + 


(1  +  6) 


Equating  the  expression  in  braces  to  zero,  and  substituting 
for  v  and  G  their  expressions  in  terms  of  x: 


v 


Si  jr 

-  x,  6  -  1  +  2Xz 

<p  m 


+  2Xx, 


we  obtain  the  possibility  of  determining  xm  for  which  the  pressure 
is  maximum: 


iu,  K  (G0  + 

8  h  sIKxm 


1  +  a 


or 


-  (i  +  e ) 


o 


i 

\ 

\ 

\  484 
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and  find  ,the  values  of  v',  /  p'.  If  p'  coincides  with  p(0\ 

it  is  indeed  the  true  maximum  pressure.  However,  if  p^  f 

then  p’  must  again  be  substituted  in  (18)  or  (19)  and  a  hew  x” 
m  m 

obtained;  then  the  whole  process  is  repeated  and  a  new  p"  is  obtained. 

m 

If  x"  is  chosen  correctly,  p"  should  not  differ  from  p’  bv  more 
m  m 

than  10-20  kg/cm^  (the  accuracy  of  a  slide  rule). 

It  must  be  remembered  that  equations  (18)  and  (19)  are  used 
for  calculating  and  can  not  be  employed  for  calculating  p^. 

When  carrying  out  the  approximations,  the  following  should  be 
kept  in  mind:  the  relation  p,  x  is  represented  by  a  curve  shown  in 
fig.  134,  which  varies  slowly  in  the  neighborhood  of  the  maximum. 

The  true  value  of  xm  is  not  known,  and  we  find  by  means  of 
equations  (18)  and  (19)  only  a  certain  approximate  value,  which, 
even  upon  substituting  the  value  p^  -  300  kg/cm2  for  in  (19) 
will  give  a  value  of  differing  from  the  real  value  by  not 
more  than  10%,  The  value  of  will  then  be  sufficiently  close 
to  the  true  value  of  pn,  and  at  the  next  approximation  x ^  will 
practically  coincide  with  x^. 

Whatever  the  quantity  of  p^  assigned  in  the  first  approximation. 

*  w 

whether  smaller  or  larger  than  the  real  value  of  pm,  the  values  of 
Pa  and  p^  will  be  smaller  in  both  cases  than  the  real  pn.  In  the 
subsequent  approximations  the  pressure  values  must  increase,  tending 
toward  the  real  pE,  l.e.,  <  P”’ —  Pn  real,  regardless  of  the 

value  of  x^0^.  Lecturer  Belenky  proved  analytically  that  in  successive 
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approximations  the  quantity  x’  is  monotonic  increasing,  tending 

toward  xm  as  the  limit,  while  x”  is  monotoaic  decreasing  and  tends 
toward  the  sane  limit  x  .  J 

to. 


Tig.  134  -  Determination  of  x^  and  pB  (graph  p,  x) . 

This  law  must  be  used  for  controlling  the  accuracy  of  the 
calculations. 

Since  the  pressure  varies  slowly  in  the  neighborhood  of 
the  maximum,  the  p^  obtained  following  the  substitution  of  these 
values  of  in  the  working  equations  will  be  very  close  to  the 
real  pm,  and  the  second  approximation  will  be  adequate  to  obtain 
a  value  of  p”  sufficiently  close  to  the  real  value. 

Having  found  xa,  we  substitute  it  into  (11)  for  v,  (10)  for 
_J14>  for  l ,  and  (15)  for  pa,  and  obtain  tne  elements  of  the 
projectile 's  motion,  i.e.,  v  u,  ^  and  pm  at  the  instant  of 
greatest  pressure. 

Expression  (18)  gives  the  analytical  expression  for  x  at 
which  the  gas  pressure  becomes  maximum.  Can  this  equation  always 
be  usjed  to  determine  the  maximum  pressure? 

In  most  cases  when  the  chosen  loading  conditions  are  normal, 

this  formula  will  give  the  right  answer.  But  there  are  cases  when 
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it  may  yield  a  value  xm  devoid  of  physical  meaning.  This  occurs  when: 


-  1  -  Zq  corresponds  to  the  instant  when  the  burning*  of  the  powder 

terminates,  the  instant  when  the  inflow  of  gases  ends.  For  this 

reason  this  formula  will  give  realistic  results  while  x  is  smaller 

m 

than  or  at  most  is  equal  to  x^(xm  <  x^.) . 

When  x  <  xw  we  have  a  normal  case:  the  maximum  pressure  is 

o  Jv 

reached  before  the  end  of  burning.  When  xE  -  x^,  the  maximum 
pressure  is  reached  at  the  end  of  burning.  Finally,  when  xn  >  x^, 
we  have  the  case  of  the  so-called  "unreal”  maximum,  i.e.,  a  purely 
analytical  case.  In  reality,  when  xm  >  x^,  the  powder,  burning 
according  to  a  definite  law,  stops  burning  on  the  upward  branch 
of  the  pressure  curve,  the  flow  of  gases  stops,  following  which 
the  pressure  begins  to  drop,  in  spite  of  the  fact  that  the  analytic 
maximum  had  not  yet  been  reached.  In  fact,  the  maximum  pressure 
in  this  case  will  be  the  pressure  p^  at  the  end  of  burning. 


Fig.  135  -  Pressure  Curve  with 
Normal  Maximum. 

A  large  value  of  xm  may  be  obtained 


Fig.  136  -  The  Maximum  Pressure 
Coincides  with  the  End  of  Burn¬ 
ing. 

when  the  parameter  of  the  loading 
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conditions,  B 
is  thin. 


S2e^/u2fu.fin  is  small; 


this  happens  when  the  powder 


GRAPHIC  NOT  REPRODUCIBLE 


Fig.  137  -  Unreal  Maximum,  x  >  x  . 

«  K 

Such  cases  occur  in  practice  when  the  firing  is  performed  with 
thin  powders  for  special  purposes. 

The  appearance  of  the  pressure  curve  in  these  three  cases  is 
shown  in  figs.  135,  136,  137. 

At  the  end  of  the  first  period,  we  have: 


"  j'.n  rx 1 

\  \  ‘  n 
/'  ■ 

ji.  , 

VI'.'*  •  Zi. 

•T’  2^7,’ 

Itoj8lh  ii 


*K  -  ,0U  - 


J. 

B, 


Xjf  •*  1  -  Zq»  lr  m  Izv.K^xK  “  ^  > 


av.  K^xK 


si 


X 


"  —  xir? Pv  "  r 


\f<*»  I  1 - Xy 

2  K 

\ 


K  K  \  8  Ux  +  /  ) 


The  same  values  will  characterize  the  start  of  the  second  period 

\ 

the  period  of  adiabatic  expansion  of  the  gases. 
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4.  SECOND  PERIOD 


The  second  period,  starting  at  the  end  of  burning  of  the  charge 
and  ending  when  the  base  of  the  projectile  passes  the  muzzle  face 
of  the  gun,  constitutes  a  process  of  adiabatic  expansion  of  the 


gases. 

This  period  is  considerably  simpler  than  the  first,  because  the 
whole  process  is  reduced  to  the  expansion  of  gases  without  the 
addition  of  energy  and  without  heat  losses. 

In  the  second  period  ^  -  1,  the  number  of  variables  is  reduced, 
the  independent  variable  is  usually  taken  to  be  the  path  l  of  the 
projectile,  and  equations  expressing  the  pressure  p  and  the  velocity 
v  as  a  function  of  /  are  derived. 

The  beginning  of  the  second  period  is  characterized  by  the 
following  data  obtained  at  the  end  of  the  first  period: 


+  -  1;  v  •  vl;  t  -  l: 


;  p  -  p, 


V  1 1 »  T-V 


The  fundamental  equation  of  the  second  period  is: 

psC^+n-fw  -  —cprnv2  -  fu>  ^1  -  , 


where 


2fu  2 

— - vn  . 

cp8m  Up 
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Since  the  gas  temperature  is  lover  in  the  second  period  than 
in  the  first,  6  should  be  made  larger  in  the  second  period,  but 
most  authors  take  an  average  value  of  G  common  to  both  periods. 

A.  Derivation  of  the  Expression  for  Pressure  in  the 
y  Second  Period 

C*  -  *x(0_7. 

The  equation  for  pressure  is  derived  from  the  adiabatic 
equation: 


w 


1+8  1+0 

pW  -  PK*K  ,  (21) 

y 

where  and  p  are  the  gas  pressures  at  the  beginning  of  the  second 

period  and  at  a  given  moment,  respectively; 

and  W  ar.e  the  free  volumes  of  the  initial  air  space  at  the 
#  • 

same  instants. 

From  equation  (21),  we  have: 


Expanding  the  quantities  and  W,  we  obtain: 

WE  “  W0  “  aw+  S^K  “  S^1  +  * 

W  "Vq  -  aw+  s(  -  s(/1  +  0. 

Substituting  these  values  in  the  equation  of  p,  we  find: 


P  -  P, 


*1  +  V 
+  / 


,i+e 
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<  0 


si 


K> 


If  we  replace  v__  by  its  expression  v„  «  - Cl  -  zrt), 

K  K  *fm  0 


2 

VK 

TT 

V 

n  p 


C2T2,..  ,  v2 

s  IK(1  -  Zq)  tfQm  _  B  e  ;i  „  \2 

— “o~ —  je;  B  t  1  ■  v  ■ 


cp  n 


and  then 


V  "  VflP 


•i + ! 


.  B8  x2 

1  -  - — (1  -  z) 

2  0 


(24) 


When  /  -  ,  we  obtain  an  expression  for  the  nuzzle  velocity  v^: 


♦ 


This  equation  is  of  great  importance  for  investigating  the  most 
desirable  solutions  when  designing  guns. 

Equations  (22)  and  (23)  or  (24)  give  the  expressions  for  the 
gas  pressure  in  the  bore  of  the  gun  and  for  the  projectile  velocity 
in  the  second  period  as  a  function  o-f  the  projectile  path  l  . 

Thus,  on  the  basis  of  the  assumptions  made,  the  equations  derived 
above  express  the  relation  between  the  conditions  of  loading  and 
the  ballistic  elements  of  a  gun  discharge  in  both  the  first  and 
the  second  periods.  They  enable  one,  for  given  loading  conditions, 
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to  compute  the  projectile  velocity  and  the  gas  pressure  at  different 
points  of  the  projectile's  motion  in  the  bore  of  the  gun,  and  to 
determine  the  maximum  pressure,  the  muzzle  pressure,  and  the  initial 
(muzzle)  velocity  of  the  projectile. 

Curves  of  p  and  v  as  a  function  of  l  will  usually  have 
the  form  shown  in  fig.-  138. 

GRAPHIC  NOT  REPRODUCIBLE 


Fig.  138  -  Normal  p,  l  and  v,  |  Curves. 

1)  Period  I;  2)  period  II. 

C.  Equations  for  Calculating  the  Temperature  of  Powder 

Gases. 

Having  solved  the  fundamental  equation  of  pyrodynamics  and 
established  the  relation  between  the  basic  elements  (p,  v,  (  ,  and  4O 
and  the  new  independent  variable  x,  and,  consequently,  also  the 
relationship  between  these  elements,  an  equation  can  be  written  for 

determining  the  temperature  of  the  powder  gases  at  any  given  instant, 

* 

t 

and,  in  particular,  at  the  instant  the  projectile  leaves  the  bore 
of  the  gun  barrel. 

\ 

The  temperature  of  the  gases  flowing  in  the  path  of  the  projectile 

♦ 

determines  whether  the  discharge  will  be  accompanied  by  a  flash, 
orjtill  be  flashless,  because  according  to  the  present  concepts  the 

flash  accompanying  a  shbt  is  a  process  involving  the  burning  of 

\ 
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inflammable  hydrogen  and  carbon  oxide  gases  making  up  about  50% 
of  the  entire  gas  mixture. 

If  the  temperature  of  these  gases  is  very  high,  the  gases 
will  burst  into  flames  when  mixed  with  the  oxygen  in  air,  and 
produce  a  flash  accompanying  the  shot. 

In  order  to  obtain  the  desired  equation,  let  us  make  use  of 
the  energy  balance  equation  in  which  Ec  is  replaced  by  R/q : 


RT«* 

~S 


spmv 


Since  RT,  -  f. 


1  - 


<fmv 


or 


1  «fne  V2  1 

2fw 


JL  Yl 


(26) 


np 


Knowing  v  and  ^  from  the  first  period,  let  us  find  T/T^  and 
then  T. 

Inasmuch  as 

slr 

v  -  — £  x,  .  . 


^  +  kjX  +  X  Xx2, 
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then,  bearing  in  mind  that: 


B  -  — 

jucf  nx 


we  get 


T  _  ^  Be  X2  1  B&  x2 _ 

T1  2  t  2  +  k1x  +  V.\x.2) 


This  equation  shows  that  the  variation  in  the  temperature  of  the 
gases  depends  upon  the  conditions  of  loading  (parameter  B)  and  upon 
the  shape  of  the  grain  (coefficients  K  and  X). 

At  the  end  of  burning  (^  "  1)  we  will  have: 

TK  '  B8  2 

JL  _  i .  _<i  .  V  - 1  -  ?v  <27> 

In  the  second  period  ^  -  1  and  we  obtain  from  equation  (26): 


T  ,  <?n0  2 

- -  «  1  -  -  y 

Tx  2 


1 


(28) 


where 


2  M  2fu> 
8p 
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At  the  instant  the  projectile  leaves  the  barrel 


ta 

~T 

vnP 


1  -  <erA. 


This  value  T^/Tj  varies  in  artillery  pieces  between  0,65  and  0.75. 
Comparing  the  value  v  .from  equation  (28)  with  the  values  T/T 

1 

aad  ^rom  equations  (27)  and  (28) ,  we  obtain  other  expressions 

for  T/T  : 


a 


( ll  + 

Vi  +  l 


e 


Tl  U  ♦  '/  ’ 


e 


yi  +  {a 


e 


This  equation  proves  that  the  temperature  of  the  gases  at  the 
Instant  the  base  of  the  projectile  passes  the  muzzle  face  depends 
on: 

1)  the  temperature  T^  of  the  burning  powder; 

2)  the  temperature  of  the  gases  at  the  end  of  burning: 


B6 

T 


(1  -  z0> 


2 


1 
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this  temperature  decreases  as  B  increases; 

3)  the  ratio  of  free  volumes  ({  +  lv)/(l,  +  [  ),  which 
depends  upon  the  path  traversed  by  the  projectile  at  the  end  of 
burning  and  decreases  as  increases. 

D.  Equations  for  Calculating  the  Time  of  Motion  of  the 

Projectile. 

The  time  t  does  not  appear  directly  in  the  solution  of  the 
fundamental  problem  of  pyrodynamics ;  one  may  compute  and  draw 
the  curves  of  the  gas  pressure  p  and  the  projectile  velocity  v 
as  a  function  of  the  projectile  path  /,  and  by  this  means  solve 
the  fundamental  problem  of  internal  ballistics,  yielding  the 
design  data  of  the  gun  (volume  of  powder  chamber,  length  of 
projectile  path). 

But  in  order  to  fully  clarify  the  phenomena  taking  place  during 
a  shot,  it  is  also  necessary  to  know  the  variation  of  the  basic 
elements  (p,  v,  /  )  as  a  function  of  the  time  t,  particularly, 
because  some  of  the  existing  devices  permit  determining  the  path  l  , 
the  velocity  v,  and  the  gas  pressure  p  as  a  function  of  the  time  t 
(velocimeter,  piezoelectric  manometer).  Moreover,  it  is  the 
pressures  curves  as  a  function  of  time  which  must  be  known  when 
solving  problems  relating  to  the  theory  of  gun  mounts  fuzes  and 
firing  devices. 

The  time  of  motion  of  the  projectile  in  the  barrel  can  be 
obtained  most  simply  if  the  curve  of  the  velocity  v  as  a  function 
of  the  path  /  is  available,  and  by  using  the  following  equation  of 
mechanics : 


61 

v  “>  — - — . 
dt 
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whence : 


If,  having  the  curve  v,  |  ,  we  plot  the  curve  -i-,  /  ,  then  by 
taking  the  integral: 

j  -H* 

o 


Ve  could  determine  the  time  of  motion  of  the  projectile  along  the 
given  path  l  .  But  inasmuch  as  at  the  lower  limit,  when  /  «■  0, 
v  -  0,  and  the  integrand  l/v  becomes  infinite  (l/v  -  1/0  «  co ) , 
it  is  impossible  to  perform  the  integration.  Therefore,  the 
time  t  is  divided  into  two  parts,  t'  and  t": 


t  -  t’  +  t". 


(29) 


where  the  first  time  interval  t’  -  from  the  start  of  motion  up  to 
a  point  representing  a  small  length  of  the  path  / -  is  calculated 
.approximately,  and  the  second  interval  t",  from  ('  to  /  along  the 
path  -  is  calculated  by  means  of  quadratic  formulas: 

t 


r 


The  first  time  interval  t’  is  found  from  the  equation: 


av. 


F-TS-7327-RE 


499 


* 


where,  in  the  first  approximation: 


0  +  v*  v, 


and  v*  is  the  velocity  of  the  projectile  at  time  t'  and  the  path 


distance  consequently: 


2  l  * 

f  -  — : 


the  smaller  the  distance  / ’,  the  greater  will  be  the  accuracy  of 
determining  t. 

Substituting  t*  and  t"  into  (105),  we  obtain  the  equation  giving 
the  time  of  motion  of  the  projectile  in  the  bore  in  the  form: 


t.Hl  + 


A.df. 

V 


Inasmuch  as  the  first  interval  of  time  for  traversing  the  path  [ 
as  determined  by  (30),  is  very  approximate,  Prof.  E.L.  Bravin 
proposed  a  more  exact  expression  for  computing  the  average  velocity 
of  the  projectile  along  the  segment  ot'.  He  assumed  the  acceleration, 
rather  than  the  velocity,  to  be  linear  along  this  segment  (fig. 

139)  : 


— <P0  +  kt), 

<j>  m  w 


where  k  »  (p'  -  p0)/t’  is  the  angular  coefficient  of  the  straight  1 i a 
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in 


P>  P'J  a  is  a  factor,  smaller  than  unity,  determined  from  the 
*  • 

condition  that  the  areas  bounded  by  the  curve  p,  t  and  the  straight 


line  P^P2  replacing  it  along  the  segment jft'  are  equal.  When 


determining  v'  /  in  terms  of  v',  the  coefficient  a  is  reduced. 


p  ’  -  p 

dv  -  (  p  + - ^ 

«pm  \  0 


&a 


t’ 


t  dt. 


/ 


Fig.  139  -  Curve  p,  t  Along  the  Initial  Path  Segment 
(According  to  Bravin). 


After  integration,  we  obtain: 


s  /  P'  -  Pn  t2 

v  «  -~—a  p  t  +  — _ 2 

«f  m  \  0  t 1  2 


Assuming  that  t  -  t\  we  find  v',  the  velocity  of  the  projectile 
at  the  time  t  *  r 


The  average  value  of  the  projectile  velocity  along  this 
•segment  is  found  from  the  following  equation: 


F-TS-7327-RE 


501 


If 


graphic  not  reproducible 


sa  1 
cpm  t* 


sa ' 
«p  in 


+  P 

- f 

6 


SQ 

or,  replacing  -  t’  from  the  preceding  equation  by  v’,  we  obtain; 

«p  m 


av. 


2P0  +  P1  v_l 
Pa  +  P*.  3 


The  final  expression  for  t*  will  be  in  the  form: 


V 


v 

av. 


3/J  Pp  +  P* 

v*  2p  +  p' 
0 


(31) 


This  is  the  equation  proposed  by  Prof.  E.,L.  Bravin/~3_7* 
Comparing  it  with  the  previous  expression  for  t',  we  note  that 
the  time  t  obtained  by  the  first  expression  is  shorter,  and  the 
difference  between  them  increases  as  the  length  of  the  segment  ( ' 
and  the  pressure  p'  increase.  At  the  limit,  when  p*  is  reduced  to 
pQ,  the  two  expressions  for  t'  become  equal: 

t,  31  •  2p0  _  21  ' 
v*  3p0  v'* 

If  the  relation  p,  t  along  the  first  segment  is  expressed  by  a 
second-degree  equation,  the  resulting  equation  will  be  more  exact: 
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,  _  4M  2Pq  +  P' 


2  v' 


5p0  +  p' 


Prof.  Bravin  also  introduced  equations  for  computing  the  time 

in  segments  measuring  0  to  t  ,  t*,  to  t„,  and  to  t  : 

®  “*  K  K  A 


3l  p_  -f  p 
t  _ jn 

n  v  „  * 

®  2 Pn  +  Pm 


3<(K  -  ,m)(pm  +  PK> 


t  -  t  «  — _ R  m  m  s 

k  "  Vp»  +  2pk)  V2pn  +  PK> 


,  *•  ‘  *k  - 


3  ((a  ~  ^kHpk  +  PA^ 

— - - - 5 - - - ~ - 

vK(pK  +  2po)  +  (2Psf  +  pn) 


A  v  r  *Vv 


In  order  to  compute  t  in  the  first  period,  the  graph  x 

P 

can  be  used  also. 

Indeed,  x  -  z  -  z^; 


whence 


*1  dt 


dt  ••  I. 
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t  - 


(32) 


The  function  under  the  integral  sign  does  not  become  infinite  when 

p0>  °- 

When  calculating  the  total  time  of  the  shot,  it  is  necessary  to 
consider  not  only  the  time  of  motion  of  the  projectile,  but  also 
the  burning  time  of  the  powder  in  the  preliminary  period  before  the 
start  of  this  motion,  which  is  computed  by  the  following  formula: 


2.303  T0 


log 


9 


1  i  A  K 

where  -  **  - -  — ,  I  «  e  /u,,  and  p  is  the  pressure  of  the 

1  -  Vs  *K  K  1  B 

igniter  gases. 

The  time  lapse  between  the  instant  the  firing  pin  strikes  the 
percussion  cap  and  the  end  of  burning  of  the  igniter  is  usually 
not  taken  into  jiccount. 


5.  SAMPLE  CALCULATION  OF  THE  GAS  PRESSURE  CURVE 
AND  OF  THE  PROJECTILE  VELOCITY  BY  THE  l  METHOD 

■  ^av. 


The  following  data  are  given: 

Barrel :  76  mm  gun,  1936  model 

Chamber  capacity,  Wq  in  dm3. . . . .,1.515 

Cross-sectional  area  of  the  bore,  including  the 

rifling  grooves,  s,  in  dm^...., . .....0.4692 

Path  traversed  by  the  projectile  in  the  bore,  }  , 

in  dm . . . . . 33.91 

Projectile : 

Weight  of  projectile,  q,  in  kg . . . 6.2 


F-TS-7327-RE 


504 


ik, -a*-**-** 


m- 


Forcing  pressure,  pQ,  in  kg/cm2 . 

*  .  „ 300 

■ 

Charge : 

Weight  of  charge,  u>  ,  in  kg . . . 

.  .  1  0ft 

* 

Powder  constants: 

Powder  energy  (force)  f,  ir.  kg*dra/kgs _ _ 

. . .950, 000 

1 

o 

Co-volume,  a,  in  drav/kg . „ . . . . .  , 

.  O  08  • 

* 

Density  of  the  powder,  &  ,  in  kg/dro3 . 

. ..1.6 

Burning  rate,  u^,  of  the  powder  when  p  -  1,  in* 
dm/sec : kg/dm2 . t 

o  fifmAfi’ ta 

Dimensions  of  strip  (thickness  2ex,  in  mm)."..... 

•  «  *  v  *  uUvUU  i  *t, 

.. .1.357 

X  -  1.06  ' 

K\-  -0.06 

Polytropic  index  k . 

1  9 

®  *  k  -  1 . . . 

0  9 

4ft 

l 

FUNDAMENTAL  EQUATIONS 

• 

• 

• 

A.  Preliminary  Period 

1  1 

A  6  / 

'V  ®*  - - :  £  ■"  \  / 1  +  i  -\tir  •  7 

T0  0  V  *  tQ *  0 

f  1  U  V  *  0  0  «(Gft  +  1) 

-  +  a  -  —  o 

P0  « 

K 

* 

B,.  First  Period 

■ 

si  tr  .  _ 

v  m  to  +  klx  +  ^AX2S 
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'■^c.-^:  *  a 


( l\  +lK 

V  ll  +  1 


v  “  vf 


T  -  Vr 


l^  m  Iq&  (  ~  ~  c  **  ~  a^)  > 


,J!l iiW.- 


fj  4-  / 


Llis' 

>x  +  '  . 


1  -  5£_<1  -  z0>2  ; 


2g  f  u> 
«f  0  q 


The  calculation  of  the  constants  is  effected  first  from  the 


given  data: 


A  «  0.7128; 


0.7128  1.6 


0.02429; 


+  a  - 


950000 

30000 


4-  0.98  -  0,625 
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4 (-0.0566)0.02423 


0.9972; 


>0  "  \/ 1  +  4  7+o 


i  + 


1.06 


2% 


0 


X(GQ  +  1) 


0.02294; 


xK  «  1  -  z  -  1  -  0.02294  -  0.97706; 


k  -  K6n  -  1.06-0.9972  -  1.057; 
1  0 


XK  "  u 


ex  0.0067S 


1  0.0000074 


916.2; 


<f>  -  1.03  +  —  —  »  1.03  +  —  .  -  1. 

3  q  3  6.2 


sl„  0.4692-916.2-98.1 

— -  -  -  -  6253; 

t?m  1.088-6.2 


1. 515 
0.4692 


3.228; 


s2l|  0. 46922*  916 . 22 • 98. 1 

B - -  - - - - -  -  2.617 

fwfB  950000*1. 08'1. 088-6. 2 


Be  2.617*0.2 

B  -  —  -KX«  — - —  +  0.06  -  0.3217 

12  2 


088; 
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COMPUTATIONS  FOR  THE  SECOND  PERIOD 
Calculation  of  the  constants  of  the  second  period: 


2  *950000 * 1 .08 *98.1 
10880.2*62 


-  149/300,000; 


-  0.250;, 

% 


?2 

C  -  1  -  -  1-0.250  -  0.750. 

,2 

nP 

Proa  the  first  period,  +  l K  “  18.384. 

-  •  • 

pK  892  kg/c»2 . 

Table  .of  the  Elements  of  the  Second  Period 


initial  Formulas 


.  i„  (tLSSi—  V  “ 

Vo.  973  +i  y 


V  -  V 


a 

Operations 

S' 

' 

Muzzle 
.  Face 

i 

1 

l 

22.60 

27.44 

‘ 

33.91 

2 

(+) 

4 

'i 

0.973 

0.973 

0.973 

3 

t+  r. 

,23.573  . 

28.413  • 

34.883 

A 

> 

• 

4 

*1  + 

u 

J-1, 

0.7799 

0.6471 

0. 5270 

'i+'  , 

1.8920 

1.8110 

1.7218 

3 

log 

-0.1080 

-0.1890 

-0.2782 

r 

-0.1296 

-0.2268 

—0.3338 

6 

|  1.2  log  | 

I. *70*  : 

T.7732 

T.6662 
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Logarithms  of  ths  Function 
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Table  (Coni’ d.) 


Initial  Formulas 

- - - « 

No.  | 

! 

Operations 

. 

4*- 

Fat  <t* 

-  12,220  y/l-0. 7501)°* 2 

7 

^1  <2 

0.7420 

0,5932 

0.4636 

♦  * 

8 

_  _  n  1.2 
p 

. 

662 

529 

414 

(■ 

9 

0.2  log  ^ 

-0,0216 

1. 9784 

r 0.0378 

1.9622 

-0.0556 

T. 9444 

10 

^0.2 

0.9515 

0.9166 

0.8798 

11 

0.750^°* 2 

0.7136 

0,6874 

0.6598 

- 

12 

1-0. 750^0*2 

0.2864 

0.3128 

0.3402 

13 

v  in  do/sec 
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6822 

1  7117 
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CHAPTER  2  -  PROP.  N.F.  DROZDOV’S  EXACT  METHOD^ 


(Written  by  Prof.  G.V.  Oppokov) 


The  assumption 


^  » 
*  Tav 


made  in  the  preceding  chapter  gives  an  approximate  solutions  Yet 
the  differential  equation  of  the  projectile  path  in  the  first  period 


dl 

dx 


Bx(/^  +  f) 


t-  'Sr* 


(33) 


can  be  integrated  exactly. 

Prof.  N.F.  Drozdov’s,  great  contribution  to  the  fieid  of  internal 
ballistics  lies  in  the  very  fact  that  he  was  the  first  to  solve  this 
equation  exactly,  without  any  additional  assumptions  or  simplifications, 
as  had  been  done  before  him  by  all  the  other  authors  without 
exception. 

Namely,  if  we  introduce  for  convenience  the  following  designation: 


M  - 


Bx 


+  -5f*2 


(34) 


equation-  (33)  takes  on  the  form: 


d|_ 

dx 


-  Ml 


-  M/^. 


(35) 
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When  this  differential  equation  of  the  ;£irst  order  is  integrated, 
the  following  relationship  obtains: 

A 


l  «■ 


f  «dx  -  ( 


Mdx 


i  e 


0  l 


Mdx. 


(36) 


/ 

len 


The  integral-exponent  of  e  in  the  .right-hand  side  of  equation 
(36)  is  (see  pp.  478-481)  equal,  as  before,  to: 


*  ,  P 

f  B-  (  ‘WP 

MX  -  -2-  _ 

J  B1  J  V  +  f>  -  p2 

n  a  • 


B 


-  InZ 


(37) 


ff : 


The  main  integral  is: 
x 

-f  Mdx 

-  X  J  X 

•  j  Mdx  -  j 


x 

f 


Mdx 


*  a.  d(-< 


). 


Noting  this  peculiarity,  the  author  integrates  by  parts: 


Y  - 


A 

-I 


Mdx 


Mdx 


0  +  J  e 

0  o 


J  e  '  .  dV 
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AWi 


* <  *  C * 


or,  ^“see* equation  (37 )_/• 


I  5 
*  ST 


T  -  -l*. 


+  •'to +  j  z .  dV 


S'  -  y0  ♦  *x*  +kXx2;  -  »t* 


vbere  a 


-  T.(-  -  4  ) 


dl^  -  -ak^dx  -  2axXxdx, 


whence 


T  -  -l^Z  *X  -  *k 


r  Bi 
o  1 


dx  -  2a*X 


B 

x  — - 

I* 


It  is  now  possible  to  proceed  in  two  ways:  eliminate,  from  the 
equation  either  the  first  or  the  second  Integral  in  the  right- 
hand  side.  The  author  selected  the  second  course. 

Nanely,  it  follows  fron  (34)  and  (37)  that: 


t  Bxdx 

l  %  +  klx  -  Bi 


r2 


In  Z 
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t 


SPi  -  ,  • 


I'ijKI.WWWIilPI/ri 


Integrating  by  parts: 


B 

B, 


xdx 


-  V(-’ -  V- - -*)  ■" 


B 


y 


B 


22 

®i 


V  >. 


B 


2B  /  B1 

-V  j  2 


xdx  + 


B 


l 


B 

®1 

Z  dx. 


because  when  x  *  0.  we  will  have : 

V  J>  *  -  V  * 


B 


B1 

1  0 

Z  -  e“  -  1. 


Te  will  find  the  desired  integral  from  equation  (38): 


B 

x  r 


-B 

B, 


f  X  2  °  1 

J  z  xdx " 772T  <B>X  ’kix'Vz 


Bi 


2£_  +  J-_  f  z1 

B  ♦  2BX  B  +  2B^  J 


dx. 


How  the  obtained  value  of  the  Integral  aust  be  substituted  into 
(38),  noting  that: 
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*,«**»*«•-  * 


V 


>«u< 


c-t/' 


I 


IH 


**>  .  i  _ 

«?T  ^  “  7tx' 


after  the  transfer  of  to  the  left-hand  side,  we  obtain: 


*  .  .  / 

..  t 


/ 


.  / 


ll  2 

If  Bi  Bi 

+  l  <b  -  c  [  Z  dx)Z  . 

o 

t 

/ 

dividing  both  sides  by  and  noting  that: 


/ 


«  .  |  _  w  _  a 


ve  obtain  Drozdov's  well-known  equation: 


(  .  /A  ¥*  *1A  /  BI  2, 

_  4  —j - j-  - - x-  I'f - — X  |  4 

0  *0  «'l  *1  V  2 


+  (b, 


•  ..  j  ■  1 


B  B 

*  r  •  §: 


dx)Z 


In  this  equation  we  have: 
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L. 


(40) 


s&m >w^w> m»ww^."«»*>- 


B 

*  r 

V  »V 


’,  £ 


(  “i  H  (  .  i 
jz  Bx-gi-  J  Z 


.The  table  for  the  last  integral  with  three  entries,  p,  ^ ,  and 
B/Br,  was  prepared  by  Prof.  D.A.  Venttsel-  and  is  reproduced  in  the 
appendix. 

Table  2  -  Computational  Formulas  Used  in  Drozdov's  Exact 

Method 


«f-  K  +  ~ 

3  q  * 


-  v 

K  u;  *'  1,0 


-  1  ; 


*±.  #  ZL. 
s»  1  s  * 


I  U  1  w  1 

A  8  &2’  *  8 


/  .1  —  m ,  L,  ^  -  f  ft  -  <  .  ‘  0  . 

tn  4  4,’  'n  s’  In 


t)=  >i- 


-  *X; 


.  Y  % 

’  •  t:4 


B  +  2B 


Clkl 


+  8,,*^) 
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(  b  Is  a  function  of  B  and  £>  for  which  a  special  table  has  been  compiled 
and  is  presented  below.  This  equation  enables  one  to  avoid 
approximations. 
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Table  3  -  Approximate  Values  -of.  the  Function 
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The  equations  of  the  preceding  table  should  be  applied  to  the 
first  period;  the  equations  for  the  preliminary  and  second  periods 
remain  unchanged  (Table  3). 

•  « 

6.  EXAMPLE  OF  CALCULATING  THE  GAS  PRESSURE  CURVE 


AND  THE  PROJECTILE  VELOCITY 
BY  PROF.  N.F,  DROZDOV'S  METHOD 

The£  following  data  are  given: 

,  Barrel :  76  mm  gun,  model  1936 
/  3 

/  Chamber  capacity,  W^,  in  dm  . . .............1.515 

Cross-sectional  area  of  the  bore*  Including  rifling, 

s,  in  dm^......  . . . . . . . . . 0.4692 

Path  traversed  by  projectile  inside  the  bore,  l  , 

in  dm . . . . . . . . f*. _ 33.91 

Projectile 

Weight  of  projectile,  q,  inkg. ....... . . . 6.2 

/  2 

Forcing  pressure,-  p0,  in  kg/cm . . . 300 

Charge 

Weight  of  charge^,  w,  in  kg . .  1.08 

Powder  constants 

Powder  energy  (force),  f,  in  kg* dm/kg. . 950000 

•  Co  volume,  a,  in  dm  /kg .  0. 98 

3 

Powder  density,  6,  in  kg/dm  . . .  1.6 

Burning  rate  of  powder,  u^,  when  p  -  1,  in 

dm/sec:  kg/dm2.. . .....0.0000074 

Dimensions  of  strip  (thickness  2e^  mm) . 1.357 


X  -  1.06 

KX  -  0.06 


Polytropic  index  k . . . . . 1.2 

e  -  k  -  l... . o.2 
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BASIC  COMPUTATIONAL  FORMULAS 
A..  Preliminary  Period 


A 


6 


4*  " 
T0 


~  +  a  T 


Po 


1  +  4 

* 


7  To’ 


zo  “ 


2*o 


*(eQ  +  l) 


B.  First  Period 


si 


-  -  x;  'V-  +  k.x  +  xXx2; 

»a  1  T0  1 


Bi 


A^+  A 


'  *i  ^  (  'v'  T*2)  +  <bi  ~  [  z  <lx>z 


P  -  fA 


,  B»  2 

T  x 


A>i>  +  A 


where 


»2I2 


BB 

B  -  - - B  - - x\; 

1  2 


o  *  *1  " 


2xX 

-  « 
2  +  2B* 
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vfikf 


X 

A 


1_ 

$ 


y 


% 


0.7128  1.6 


i  +  awi 

p0  a 


950000 

30000 


-  0.02429; 


+  0. 98^0. 625 


s 

1  • 


6o  " 


.  r  4 (-0.0566) 0.02429 

1  +  4  — -  \/l  +  - — - —  -  0.9972 

K  ^0 


1.06 


Z0  ’ 


2+o 


2-0.02429 


*(*0  1)  1.06*  1.9972 


-  0.02294; 


Xg  •  1  -  z*  «■  1-0.02294  ■  0.97706; 


kl  "  **6o  “  1*06*0.9972  -  1.057; 


e.  0.00678 

- -  -  916.2; 


K  u 


1  0.0000074 


1  w  1  1 • 08 

-  1.03  * - -  1.03  +  —  -rr—  -  1.088; 

3  <!  3  0.2 


l 


si *  0.4692*916.2*98.1 

—A  -  -  -  6253; 

1.088*6.2 


1.515 

0.4692 


3.228 
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c.k.  0.2576  -  1..  057 

JLi.  -  - - - -  -  0.8464; 

B1  0.3217 


0.3217*0.02429 

1.0572 


0.006995; 


2*617 

- — —  -  8.134; 
0.3217 


f-A  -  950,000-0.7128  -  677,200. 


Computing  xa 


x  - 


B  +  2B, 


+  b. 


<♦: 


-a 


a 


where 


h 
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e  for  Computing  by  Means 
Equation  (*) 
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Table  for  Computing  the  Ballistic  Elements  for  the  First  Period 

(Calculated  on  a  50-cm  Slide-Rule) 


006995  15 


This  denotes  the  page  number  in  the  original  manuscript.  This  portion  has  not  been  translated  as  yet. 
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COMPUTING  THE  SECOND  PERIOD 
Conputation  of  the  constants  of  the  second  period: 


From  the  first  period  +A^  *  6.1028,  -  832.6  kg/cm 


F-TS-7327-RE 


534 


for  horizontal  interpolation; 


jfn+1  ~  U 


- 


P  -  B 


'ft+1  “  pn 


for  vertical  interpolation. 


Example.  Find  log  Z”1,  if  it  is  known  that  J*  -  0.006995;  p  -  0. 
These  values  of  y  and  P  are  not  found  in  the  table  of  the  logarithms 
of  function  Z  *,  and  we  take  the  nearest  values,  i.e. : 


f  “  0.006  and  ^  -  0.008; 


P  «*  0;28  and  p  -  0.300, 


and  we  find  the  values  of  log  Z"1  for  these  combinations.  We  then 
calculate  and  Cp 


0.4975 

r 


P  X  0.006  0.006995  0.008 

0.28  0.1309  0.1295  0.1280 


0.2973 


0.1397 


0.300  0.1428  0.1413  0.1398 


0.006995-0.006 


0.008-0.006 

0.4975 

0.2973-0.28 


0.300-0.28 
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We  find : 


log  Z“1}  -  0.1309  -  0.4975(0.1309  -  0.1280)  -  0.1295; 


log  -  0.1428  -  0.4975(0.1428  -  0.1398)  -  0.1413. 


_Fiaaily  »e  obtain: 


log  Z”1  -  0.1295  +  0.865  (0.1413  ~  0.1295)  -  0.1397. 
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DETERMINATION  OF 


f  B/Bj 

Iz  ' 


dS  FROM  THF  TABLES 


r  ,B/Bi 

r  ‘ 


(APPENDIX  3) 

Prof.  N.F.  Drozdov  recommends  to  calculate  the  quantity 


by  the  trapezoid  rule.  •  In  order  to  simplify  the  'calculations,  \Z  dx 


may  be  found  from  the  tables  of  the  function: 


D  £ 
,  B, 


where 


whence 


\p  -  r~*> 
\  i 


Consequently 


x  1  p  JL\ 

C  *l  ki  Y  Bi\ 

|  Z  dx  -  g-~  Z  dD. 


P.  ~\ 


Z  dp. 


—  \  f  B/®1  .  / 

The  tables  of  the  function  Z  d£  are  computed  for  B/B.  - 

\  J  1 
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-  5,  6*  7,  8,  9  and  10. 

In  our  problem  the  ratio  B/B^  -  8.134,  i.e.,  it  is  intermediate 

■'  ir" 

between  B/B,  »  8.0  and  B/B.  -  9.0.  This  requires  an  additional  * 
1  1  p 


interpolation.  Therefore  the  determinatio 


n  of  j 


Z 


dp  when  B/B^  - 


■P  ,  0 

(  B/Bl 

■  8.134  is  reduced  to  the  calculation  of  1  Z  d£  from  the  tables, 

0 

first  when  B/B,  «•  8.0,  then  when  B/B-,  -  9.0.  Interpolating  these 

,  £  B/Bi 

values  of  the  integrals  along  B/B^,  we  finally  obtain  j  Z  Ad?  for 

0 


B/b^  «  8.134  and  the  given  values  of  fi. 

It  should  be  remembered  when  performing  these  calculations 
that  the  interpolation  of  the  intermediate  integrals  must  be 
performed  by  the  same  procedure  *s  that  of  the  function  log  Z"*1. 
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0.0t>6  0.006995 


FORM  FOR  DETERMINING  \  Z  FROM  THE  TABLES 


0.085  '  K).I20  0.080 


0.0845  0.085  ^  [0.120  0.080  0.0815  O.Q83 
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CHAPTER  3  -  SOLUTION  OF  THE  PROBLEMS  OF  INTERIOR  BALLISTICS 

FOR  THE  SIMPLEST  CASES 

1.  SOLUTION  OF  THE  PROBLEM  FOR  THE  CASE  OF  INSTANTANEOUS 

BURNING  OF  THE  POWDER 

the  colloidal  powders  now  used  burn  gradually,  in  parallel 
layers,  and  when  the  web  thickness  is  properly  selected,  permit 
the  regulation  of  the  flow  of  gases  during  burning,  so  that  the 
maximum  pressure  in  the  bore  pB  would  not  exceed  a  given  value 
(usually  of  the  order  of  2500-3500  kg/ca2). 

The  case  of  the  instantaneous  burning  of  the  charge. is 
anamalous  and  generally  does  not  occur  in  practice.  It  can  be 
achieved  in  practice  only  under  special  conditions,  such  as,  for 
example,  when  burning  a  charge  of  dry  pyroxylin  in  powder  form, 
or  of  fine  porous  powder  loaded  very  densely. 

In  that  case,  if  the  loading  density  were  normal  (A  -  0.50-0.75), 
the  pressure  prior  to  the  projectile's  displacement  would  reach  a 
maximum  value  of  the  order  of  several  tens  of  thousands  of  atmospheres 
(20,000  to  40,000  kg/cm2).  The  present  ultimate  strength  of  gun 
barrels  is  such  that  the  walls  of  the  barrel  would  burst  when  subjected 
to  such  pressures. 

Nevertheless,  the  case  of  the  Instantaneous  burning  of  a  charge 
is  very  interesting;  its  examination  has  an  important  meaning  when 
compared  with  gradual  burning  of  powder  because  in  so  doing  the 
Importance  of  slow  burning  and  of  the  shape  and  dimensions  of  the 
powder  grains  become  evident.  Moreover,  the  pressure  curve  p,  / 
in  the  case  of  instantaneous  burning  becomes  a  sort  of  a  "guide" 
for  the  curves  depicting  slow  burning.  These  p,  /  curves  arrange 
themselves  with  a  certain  regularity  with  respect  to  the  instantaneous 
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burning  curves. 

The  analytical  solution  of  the  problem  is  very  simple  in  the 
case  of  instantaneous  burning,  because  one  of  the  four  variables 
.  entering  into  the  fundamental  equation  of  pyrodynaraics  is  transformed 
^  into  a  constant  ('V  -  1). 

Let  the  gun  and  the  loading  conditions  be  characterized  as 
lollows: 

The  chamber  capacity  is  the  cross  sectional  area  of  the 

bore,  including  the  rifling  is  s,  the  path  of  the  projectile  is  t  , 

the  weight  of  the  charge  iso,  and  the  weight  of  the  projectile  is  q. 

The  energy  of  the  powder  is  f,  and  a  is  the  covolume;  the  adiabaitic 

index  is  k  -  1  +8,  and  the  secondary  work  done  is  taken  into- 

account  by  the  coefficient  tp «  a  +  b  — . 

Q 

When  ~  1,  the  fundamental  equation  of  pyrodynaoics  is: 

ps(  +  l)  ■  fw-  -~-<pmv2;  *  (41) 

the  equation  of  the  projectile  motion  is: 

psdl  -  mvdv,  (42) 

* 

where  /  ■  (W0  -  au)/s  is  the  reduced  length  of  the  chamber  at  the 

end  of  burning. 

When  the  powder  in  the  chamber  is  burned  instantaneously  the 
maximum  pressure  is  determined  by  means  of  the  well  known  formula: 
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(43) 


✓  - - - - : 

1  1  -  at,  rQ  -  aw  si 


The  projectile  will  be  set  in  motion  when  the  following  initial 
conditions  obtain: 

/ 

ft 

i 

l-  0;  v  -  0;  p  -  PV  •  / 

Eliminating  pressure  p  from  equations  (41)  and  (42),  we  obtain 
v  it  i  function  of  / : 


<Jl  '  q>mvdv 

li  4-  f  a  » 

1  fu  -  — 

2 


t«fmv^ 


-  0  2 
fw  -  -^mv 


We  shall  integrate  this  differential  equation  with  the  variables 


■eparab'  .  ; 


l\+l 


0<pmv^\ 


-  1,1  - 


(44) 


whence 


‘  ’n) 


(1  +  ( 
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/.  -  ■>  ■*  •*  -*WV»  -- 


-  - 


itk 


whore 


2 


v 


np 


2gfu»/«pGq  is  the  limiting  velocity  of  the  projectile: 


(45) 


4 


This  formula  expresses  the  velocity  v  of  the  projectile  as  a 
function  of  its  path  t  ;  the  velocity  increases  when  /  increases. 

In  order  to  determine  the  dependence  of  p  on  /  from  (44)  we 
determine : 


fw 


6  2 

— <|>mv 
2 


fw 


u  +i 


i 

» 

r 


and  include  it  into  (41) : 


ps(/1  +  /  ) 


Whence 


P  - 


t  U 


1+8 


■7i  V'i 


+ 1 


-  Pi 


( 


But 


wo  -  a~ 


/1  +  /  W 1  +  Bi  -  OU  +  S  l 


(46) 
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is  the  ratio  between  the  free  volumes  in  the  initial  air  space 
■eisured  at  the  start  of  motion  and  at  the  given  instant.  Consequently, 
formula  (46)  is  the  equation  of  the  adiabntic  curve  starting  with 
the  Motion  of  the  projectile  under  the  pressure  p^  -  fA/(l  -  aA) . 

The  change  in  temperature  of  the  gases  doing  the  work  is 
expressed  by  the  following  relationship  for  the  adiabatic  process: 


If  we  divide  the  numerator  and  denominator  in  parentheses  in. 

\ 

formulas  (45)  and  (47)  by  and  designate  l/l^  by  y,  we  will  get: 


(48) 


The  quantity  y  is  the  ratio  of  the  relative  projectile  path 

\ 

to  the  reduced  length  of  the  free  volume  in  the  chamber  at  the  end 

\ 

of  burning,  and  is  chlled  the  "number  of  free  volumes  of  gas 

—  \ 

expansion."  \ 
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Equations*  (48)  and  (49)  show  that  under  the  given  loading 
conditions  (q,,uj,  f ,  a,  W0,  s)  the  pressure  p  and  the  velocity  v 
depend  only  upon  the  number  y  of  free  volumes  of  expansion.  The 
greater  y,  the  higher  is  the  projectile  velocity  and  the  smaller 
the  pressure;  the  greater  the  reduced  length  1^  of  the  free  space 
in  the  chamber,  the  greater  will  be  the  gas  pressure  for  a  given 
projectile  path.  Consequently,  the  drop  in  pressure  as  a  function 
of  the  projectile  path  will  be  slower  in  a  large  chamber  than  in 
a  small  one. 

It  can  be  proved  that  the  velocity  of  the  projectile  computed 
by  means  of  formula  (48)  for  the  case  of  instantaneous  burning 
will  be  always  greater  than  the  true  velocity  for  the  case  of  slow 
burning,  under  the  same  charging  conditions. 

Indeed,  the  maximum  work  done  by  a  powder  charge  ^  of  energy 
f  in  setting  a  projectile  of  mass  m  in  motion,  is  determined  by 
the  expression  fw/8 .  This  maximum  work  will  be  the  same  for  both 
modes  of  burning  (instantaneous  and  gradual)  and  is  expressed  by 
the  areas  under  the  curves  sp  as  a  function  of  (  ,  when  /  varies 
between  0  and  infinity.  Consequently,  in  both  cases  the  areas  will 


be  equal  to: 


pd!  “ 


In  the  case  of  instantaneous  burning  the  curve  p,  [  starts 
from  the  maximum  pressure  p^,  then  varies  according  to  the 
adiabatic  law,  decreasing  continuously  (fig.  141,  curve  I),  When 
burning  is  gradual,  the  curve  II  of  the  pressure  p,  (  rises 
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gradually  from  p^,  losing  a  portion  of  area  A;  and  inasmuch  as 
the  total  area  under  the  curve  p,  /  in  the  second  case,  limited 
by  (  -  od  f  must  be  the  saue  as  the  first,  curve  II  must  necessarily 
cross  curve  I  during  burning  of  the  powder  when  the  pressure  drops, 
and  then  continues  to  rise.  The  excess  area  B  between  the  curves, 
when  t  -  oo is  the  limit,  must  be  equal  to  A.  But  inasmuch  as  the 
actual  bore  has  a  finite  projectile  path  /^  ,  the  portion  of  the  area 
B  on  this  finite  length  is  always  smaller  than  A,  and  consequently, 
for  a  given  path  length,  the  work  done  by  the  gases  and  the  velocity 
of  the  projectile  will  be  always  smaller  in  the  case  of  gradual 
burning  than  in  the  case  of  instantaneous  burning. 


Fig.  141  -  Curves  p, t  and  v, (  Depicting  the 
Instantaneous  Burning  of  Fowder. 

The  actual  initial  (muzzle)  velocity  of  a  projectile  of  a 
’medium-caliber  gun,  measured  experimentally,  represents  80-90% 
of  the  velocity  computed  by  formula  (45). 

Inasmuch  as  the  work  represented  by  the  area  under  the  curve 
p,  I  is  larger  in  the  case  Of  instantaneous  burning  than  in  gradual 
burning,  especially  at  the  Start  of  the  motion,  the  corresponding 
velocity  curve  rises  more  steeply  at  first.  Thereafter,  because  of 
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graphic  not  reproducible 


the  addition  of  the  area  B,  the  velocity  increase  becomes  greater 

in  the  case  of  gradual  burning,  and  curve  II  gradually  begins  to 

/  , 

approach  curve  I  (see  fig.  141),  tending  at  /  «  co  to  the  common 
limit  Yflp  - 

2.  SOLUTION  OF  THE  PROBLEM  FOR  A  POWDER  WITH  A  CONSTANT 
BURNING  AREA,  WHEN  THE  FORCING  PRESSURE  IS  ABSENT. 

The  solution  of  the  problem  of  internal  ballistics  for 
degressive  powders  in  the  presence  of  forcing  pr/ssure  results' 
in  equations  which  do  not  give  an  immediate  relation  between  v, 
p,  and  /,  and, " therefore, exclude  the  possibility  of  an  analytical 
examination  of  the  basic  relations.  In  order  to  obtain  this 
possibility,  it  is  necessary  to  introduce  certain  simplifications 
into  the  initial  data,  namely: 

1)  Consider  a  powder  having  a  constant  burning  area 

*  *•  1;  A-  0;  »  z. 

2)  Consider  the  forcing  pressure  to  be  negligible;  assume 
that  the  projectile  is  set  in  motion  when  the  pressure  equals  the 
pressure  of  the  igniter  gases,  and  that  the  burning  of  the  charge 
begins  when  the  projectile  is  set  in  motion: 

•>0  ■  Pb!  to  -  «• 

3)  Assume  that  a  »  l/<5 . 

When  these  assumptions  are  made,  the  solution  of  the  fundamental 
system  of  equations  is  greatly  simplified. 

The  first  assumption  corresponds  to  the  burning  of  long  tubular 
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powder;  the  second  corresponds  to  projectiles  with  pre-cut  bands; 
the  third  assumption  simplifies  the  solution  and  permits  determining 
the  qualitative  effect  of  the  loading  conditions. 

Under  the  assumptions  made,  the  preliminary  period  does  not 
exist*  The  motion  of  the  projectile  begins  under  the  following 
conditions: 

-  p_ ;  y  »  0 ;  z  «  0;  l  -  0 ;  v  -  0 . 

*  0  B  'o  0 

Inasmuch  as  a  -  1/8 , 

(.-'t  '  ?1  ■  ?0a  '  ai)- 

-  The  law  governing  burning  of  powder,  ^  -  f(z)  will  be 
expressed  by  the  formula: 

t  -  z  -  x,  (50) 


and  may  be  taken  as  the  independent  variable.  Then  the  equation 
of  the  projectile  velocity  will  take  on  the  form: 


v  - 


si 

J 


K  , 

it 


(51) 


The  fundamental  equation  of  pyrodynamics  is; 


ps(/^  +/  )  **  fu»*f  -  -^-cpmv^  *  fu> 


The  equation  of  the  elementary  work  done  is: 
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psdi  -  q>mvdv 


<fW. 


't’cH; 


K 


the  differential  equation  of  the  projectile  path  will  be: 


_dj _ 

(l  +  / 


Bd^ 


£ 

0 


Integrating,  we  get: 


1 


whence. 


(52) 


l  - 


(53) 


Designating,  as  in  the  case  of  instantaneous  burning. 


we  obtain  from  equation  (52) : 


t 


\ 


1 


1 


T 


(1  +  y)  2 
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(54) 


\ 


v 


<&»  -or  «\  .  .  5? 


and  inasmuch  as 


V  ~  - vl/. 

n  1 


Sl 

(l  +  y)  2 


(55) 


These  equations  give  the  direct  dependence  of  v  and  on  the 
path  of  the  projectile  l,  or  y  « 

Let  us  write  the  pressure  equation,  taking  into  account  the 
igniter  pressure: 


fBwB  fBAB 


p  -  = — -  -  - ; 

wo-f  wo~°~  8*1 


Vb  +  -  Tfmy2 

8(^+0 


Let  us  designate  the  relative  energy  of  the  igniter  gases  by: 


-v 

ab* 


Carrying  outside  the  parentheses  in  (56)  and  replacing  v  according 
to  (52),  we  obtain: 
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XB  +)']!  -  —  * 

si,  (1  +  y) 


1  +  y 


+  PI 


(l  +  y) 


where  «  fu>/s/^  «•  fA/(l  -  cxh)  is  the  maximum  pressure  developed 
by  the  burning  of  the  entire  eharge  within  the  space  of  the  chamber 
when  the  density  of  the  loading  is  A ; 


Be  i 

1  -  - — - - 2T  on  the  basis  of  equation  (52). 

2  JL 


(l  +  y)  2 


Substituting  -P  by  its  expression  in  (54),  we  obtain  the  pressure 
p  as  a  function  of  the  path  of  the  projectile: 


P  -  P 


a8  ,  2 

!— +Pl*nr  1 


(1  +  y)  2  (l  +  y) 


1+f 


The  quantity 


*B  fBu,B  _  PB,0 

1  +  y  81^(1  +  y)  l  +  y 


represents  the  pressure  developed  by  the  igniter  gases  in  the 
variable  space  of  the  bore.  At  the  start  of  motion  y  -  0,  p„  •*  pD  n; 

*8  the  projectile  moves  forward  and  y  increases,  pg  decreases  (1  +  y) 

times,  and  may  be  neglected  when  compared  with  the  pressure  developed  by  the 
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V 


gases  of  the  powder  charge. 

^  ' 

At  the  beginning  of  motion,  y  -  0,  the  second  term  is  equal  to 

zero  (equation  58) ;  P  “  Pg  q-  As  the  p/ojectile  moves  and  y  increases, 

*  -(i+e/2) 

.the  factor  in  the  brackets  increases,  while  the  factor  (1  +  y) 

decreases. 

The  maximum  pressure  pR  will  occur  at  some  value *y  . 

Let  us  designate:  ' 


F(y)  - 


(1  +  J)  2  (1  +  y) 


(l  +  y) 


-U+e) 


Differentiating  F(y)  with  respect  to  y,  and  equating  the 
derivative  to  zero,  we  find: 


1  +  y  - 

JM 


i  +  e 


"  ^^(8)  -  const. 


/"when  8  «  0.2;  Fi<0)  -  2.387  7 


whence. 


-  (1  C*!®)  -1_7  -  ioa  -  a  A  )/~F1(8)  -  1J 


Substituting  (59)  in  (58),  (54)  and  (55),  *e  obtain  the  expressions 
for  all  the  elements  of  motion  at  the  instant  p  -  d  : 


max. 
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fui  '  1 

si  l  +  e  * 

K 


(63 


Equations  (60^  (61),  (62)  and  (63)  give  the  direct  relationship 
between  the  elements  of  motion  and  several  characteristics  and 
parameters  at  the  instant  of  maximum  pressure.  Thus  the  path  l 

BS 

is  directly  proportional  to  the  reduced  length  of  the  chamber  (q 
and  to  1  -  aA.  When  A  increases,  lm  decreases.  The  pressure  pa 
is  directly  proportional  to  the  pressure  p^  of  instantaneous  burning 
determined- by  Nobel’s  equation;  pw  is  inversely  proportional  to 
Prof.  Drozdov's  parameter  B.  is  also  inversely  proportional  to 
the  parameter  B. 

When  ©  -  0.2: 


1.387(1  -  aA). 
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The  equations  are  very  simple  and  accessible  to  analysis. 
At  the  end  of  burning*  ^-1. 

sl» 

vr  -  — ~ ; 

A  «p» 


(64) 


i  +  jt- 


(■-“) 


J. 

e 


(65) 


(l  +  yr) 


p,  1 

1+  -I-  1 


B9 

2 


2+0 

"e~ 


(66) 


’^n  the  second  period  we  find  the  following  relationships: 

I 


*1  +  *K 

P  “  Pr  1 - 

'  K  w,  +  / 


1+0 


-  Pi 


.i+e 

\i  +  y  1 


Substituting  for  pK  and  1  +  their  expressions  from  (65)  and 


(66),  wd  find: 


P  “  - 


1  - 


B0 


(1  +  y) 1+0 


(67) 


This  i ts  the  equation  of  an  adiabatic  curve  with  initial  ordinate 

BJL 
2  r 


Pj/<1 - the  latter  increases  with  B,  i.e.,  with  the  thickness 
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of 'the  powder. 


From  the  general  equation  for  the  projectile  velocity,  we  have 


The  Temperature  of  Powder  Gases. 
In  the  case  of  instantaneous  burning  we  had: 


1 

(l  +  y)*’ 


When  burning  is  gradual: 


(61 


T 

T1 


1  _  _L 
* 


1  -  ^  v 


nP 


t  v?  • 

■  nP 


Substituting  the  values  of  v  from  (51)  and  (55)  and 
2fw/^0m,  we  get: 


T 

Ti 


1  - 


sIK  ?6m2 iw 

2fwesi 


K. 


1  - 


(l  +  y)  2 
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fK«»-  ^ .-y~  k 


(i  +  y)  2 


A  comparison  of  (69)  with  (70)  will  show  that  in  the  expression 
for  relative  gas  temperature  for  gradual  burning  of  a  powder  with 
a  constant  area,  the  value  of  the  exponent  is  one-half  of  that 
of  instantaneous  burning.  Consequently,  the  temperature  drop  in  the 
case  of  gradual  burning  proceeds  almost  at  half  the  rate  of 
Instantaneous  burning,  because  the  work  done  in  traversing  a  given 
path  is  considerably  smaller. 

The  expression  for  the  temperature  T/T^  may  be  written  as  a 
function  of  ^  only? 

.2*2  v2 

_I_  .  i  .  JV  .  i  _  _5_t  -  i  .  S».+,  (7v, 

Hr  %  *nP  2 


(71) 


i.e.,  the  temperature  of  the  gases  inside  the  barrel  during  burning 
of  powder  with  constant  area  is  a  linear  decreasing  function  of 

At  the  end  of  burning  (f  ~  1)  • _ 


-  1  - 


The  thicker  the  powder,  the  larger  is.B,  and  the  lower  is  the 
temperature  T^. 

In  the  second  period,  from  expression  (68) : 
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-A 


i  1 

!-»(!+  y,9 


(72) 


or 


(73) 


These  relatively  simple  equations  enable  one  to  perform  an 

/ 

analysis  of  the  variation  of  the  elements  of  a  shot  (p,  v,  ^  ,  T) 
as  a  function  of  y  -  the  relative  path  of  the  projectile  ~  and  to 
establish  a  series  of  relations  and  properties' of  the  variation 
curves  of  these  elements. 
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CHAPTER  4  -  ANALYSIS  OF  TS1E  BASIC  RELATIONS  FOR  THE 

SIMPLEST  CASE 

(*  •  *’  ’O  -  °>  «  -  T  ) 

* 

1,  ANALYSIS  OF  THE  FUNDAMENTAL  CURVES  p>  v,  T-,  ^  . 

An  analysis  of  the  equations  obtained, and  of  the  curves  represented 
by  them,  shows  that  they  represent  certain  staple  combinations  of  two 
types  of  curves  (fig.  142) : 

a)  Two  polytropic  curves  with  exponents  k  «  1  +  9  and  k*  -  1  +0/2, 
starting  from  the  point  (/l,  0),  first  dropping  steeply  and  then  more 
gradually; 

b)  Two  curves  analogous  to  the  polytropic  curves  but  with 
exponents  smaller  than  unity  (k  -  1  «  0  and  k'  -  1  «  0/2)  also  starting 
from  point  (1,  0)  and  descending  much  more  slowl>  with  the  convex 
side  directed  downward. 


Fig.  142  -  Basic  Types  of  Curves  for  the  Simplest  Case. 
Indeed,  in  the  case  of  instantaneous  burning: 
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!l„_j _ 

pi  i+^ 

(i  +  yK>  2 


In  the  second  period; 

_ _ 1 

Hi  +  y>1+'" 


(66’) 


(67’) 


(72) 


(68’) 


Etch  of  these  equations  contains  one  of  the  polytropics  indicated 

above  in  the  fora  of  a  variable  component, 

T 

The  curves  — =~,  y  for  instantaneous  and  gradual  burning  are 

Ti 

■expressed  directly  in  the  first  period  by  curves  with  exponents  6 
and  ®/2,  according  to  equations  (69)  and  (70). 

« 

The  ordinates  of  the  curves  ^ (54  * )  and  v  (55’)  are  obtained 
1  from  the  ordinates  of  the  curves  A T/T^  -  1  -  T/T^  (fig.  142),  measured 

from  the  horizontal  1-1,  multiplied,  by  the  coefficients  2/B9  and  2fo/lK0s, 
respectively,  and  laid  off  upwards  along  the  abscissa  from  the  origin 
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of  coordinates  to  y^  which  corresponds*  to  the  end  of  burning  (fig. 

V 

143).  The  curves  y,  y  and  - ,  y  are  inverted  with  respect  to  the 

AT  vnp 

curves  ,  y. 

1 

The  ordinates  of  the  relative  pressure  curve  (57')  are  obtained 

by  multiplying  the  ordinates  of  the  ^  curve  (54’)  by  those  of  the 

1+9/2 

auxiliary  polytropic  curve  y  -  l/(l  +  y)  .  This  same  polytropic 
curve  is  the  geometrical  locus  of  the  pressure  Pj^/pj  at  the  end  of 
burning,  expressed  as  a  function  of  y  (66). 


Fig..  143  -  Curves  ^  and  v/v  in  the  First  Period. 

It  is  seen  from  equation  (66)  that  when  B  (powder  thickness) 
increases,  decreases,  while  yK  increases  /“according  to  65)_/. 

The  curves  of  the  second  period  elements  start  to  the  right 
of  the  abscissa  y^;  the  curve  T/T^  represents  the  curve  T/T^  -  1/(1  +  y) 
of  instantaneous  burning,  multiplied  by  the  quantity  1/(1  -  B0/2) > 1, 
which  means  that  the  gradual  burning  curve  is  higher  than  the 
instantaneous  burning  curve,  the  difference  in  height  increasing 
with  the  parameter  B,  i.e,,  the  difference  being  greater  for  thicker 
powders. 

From  equation  (60)  it  is  seen  that  the  maximum  pressure  is 
independent  of  the  parameter  B,  of  the  thickness  of  the  powder,  and 
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of  the  weight  of  the  projectile.  For  a  given  A,  the  length;  l  is 

m 

proportional  to  the  length  7^  of  the  chamber;  when  A  increases  the 
imaximum  is  shifted  toward  j/ie  start  of  motion. 

/  The  real  maximum  pressure  /“equation  (61)^7  is  proportional 
to  the  maximum  pressure  in  the  case  of  instantaneous  burning, 
p^  **  fd/(l  -  q&),  and  is  inversely  proportional  to  the  parameter  B, 

or  to  the  powder  density  squared.  When  'the  energy  of  the  powder, 

;  .  /  ' 

entering  in  the  expressions,  for  the  pressure  p,  and  the  denominator 
B,  is  changed,  the  maximum  pressure  varies,  proportionally  to  the 
energy  of  the  powder  squared,  and  to  the  weight  of  the  projectile. 


because 


.2i£g 

fw<fq 


T  r2(9) 


fA  f<-T* 

F  (9)  -2~2. 

2  1  -  an  szl£ 


When  B  increases,  and  vM  decrease  also  /“equations  (62)  and  (63)^7* 

The  adiabatic  pressure  curve  in  the  case  of  instantaneous  burning, 
1+9 

p/p^  -  1/(1  +  y)  »  acts  as  "guide" . for  the  adiabatic  curves  of 
the  second  period  when  the  powder  burns  gradually.  The  ordinates 
of  these  curves  are  obtained  by  multiplying  the  ordinates  of  the 
first  curve  by  l/(.l  -  B8/2)  >  1  /“equation  (6T)'/7*  The  thicker 
the  powder,  the  larger  is  B,  the  larger  is  this  value,  and  the 
higher  will  the  adiabatic  curve  of  the  second  period  lie  above  the 
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adiabatic  curve  of  ins-tantaneoue  burning.  The  ratio  of  the  ordinates 
of  these  adiabatics  for  the  same  value  of  y  (or  I  )  is  constant 
and  equal  tp  1/(1  -  B&/2)  -  const. 

In  the  case  of  instantaneous  burning  and  in  the  second  period, 
the  projectile  velocity  is  proportional  to  the  square  root  of  the 
temperature  drop,  rather  than  to  the  first  power  of  this  factor, 
as  is  the  case  for  the  first  period. 

For  a  given  charging  density  and  during  the  burning  of  the 
powder,  the  projectile  velocity  in  a  given  section  does  not  depend 
upon  the  weight  of  the  projectile.  Indeed,  from  equation  (55) 


it  is  seen  that  for  a  given  value  of  y  -  l/lQ(l  -  oa)  and  for  one 
and  the  same  powder  (f,  a,  1^)  the  projectile  velocity  is  independent 
of  the  weight  q  of  the  projectile.  The  same  may  be  said  of  the 
temperature  of  the  gases  /"according  to  formula  (70)_7» 

If,  all  other  conditions  being  equal,  we  vary  only  q  which 
enters  into  parameter  B,  the  pressure  p  and  ^  from  equations  (54) 
and  (56)  and  the  maximum  pressure  pB  and  increases  in  proportion 

with  q,  while  the  location  of  the  maximum  and  of  the  value  v  does  not 

is 

change. 

Consequently,  the  velocity  curves  coincide  point  by  point 

when  superimposed  on  each  other (*),  and  only  when  the  projectile 

(*)  The  result  obtained  (when  a  -  1/5)  can  be  confirmed  by  comparing  it 
with  the  GAU  tables  (ordnance  tables),  compiled  for  pQ  -  300,  a  -  1,  and 
*  «  1.06. 
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Is  heavier  will- the  velocity  v  /“from  equation  (64)_7  be  attained 
earlier. 


Fig.  144  -  Ballistic  Curves  for  the  Simplest  Case. 

The  fact  that  the  velocity  of  the  projectile  does  not  vary 
with  its  weight  may  be  used  to  verify  the  complete  burning  of  the 
powder  in  the  gun.  If  a  gun  using  the  same  charge  and  type  of 
powder  is  used  to  fire  two  projectiles  of  different  weights  and 
the  velocity  remains  unchanged,  it  is  proof  that  the  combustion 
of  the  powder  was  incomplete  in  both  cases. 

Figure  144  illustrates  the  basic  curves  of  the  elements  of 
a  shot  ty,  v,  T,  p)-  as  a  function-  of  l  or  y  for  the  cases  of 
instantaneous  and  gradual  burning  of  the  powder. 

All  the  curves  on  the  graph  are  marked  with  the  number  of 
the  equation  they  represent. 

There  are  .two  basic  points  on  the  ordinate  axis,  one  is  at  - 
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and  the  other  at-l/(l  -  B8/2). 

The  two  pressure  poly tropic  curves  with  exponents  k  and  k', 
and  the  two  "polytropics"  of  temperature,  with  exponents  8  for 
instantaneous  burning  (69),  and  0/2  for  gradual  burning  (?0)  issue  from 
the  first  point;  the  ordinate  at  y  r  ~1  is  =a  common  asymptote  for 
all  of  these  curves. 

The  curves  p/p^  (67V)  and  T/T^  (72)  for  the  second-period, 
issue-  from  the  second  point,  whose  ordinate  is  1/(1  -  Bw/2).  These 
curves  are  real  only  to  the  right  of  the  ordinates  p^/p^  and  T  /T^  -  1  - 
-  B6/2  with  abscissa  y^. 

Both  of  these  curves  lie  above  the  corresponding  curves  for 
instantaneous  -burning,  the  ratio  between  the  two  sets  being  constant 
and  equal  to  l/(l  -  B8/2)  >1.  The  horizontal  line  whose  ordinate 
equals  unity  is  the  origin  for  the  curves  p/p^l)  and  T/Tj  (69) 
and  (70),  the  terminal  point  for  and  v/vK  (54*),  and  is  an 
asymptote  for  v/v^  for  both  instantaneous  and  gradual  burning  of 
powder. 

The  curves  ,  y  and  y  in  the  first  period  are  similar  to 

vnP 

the  curve  AT/t- which  is  measured  from  the- horizontal  along  the- 
oiviisnate  equal  to  unity. 

The  curve  p/p^  for  gradual  burning  is  obtained  by  multiplying 
the  ordinates  of  curve  vy(54’)  and  of  the  adiabatic  curve  2  with  the 
exponent  k*  •  1  +  8/2. 

2.  THE  CONDITIONS  FOR  MAINTAINING  THE  MAXIMUM  PRESSURE  CONSTANT. 

This  question  is  very  important  in  the  ballistic  design  of 
guns,  because  the  condition,  generally  imposed  is  that  pB  must  not 

* 

exceed  a  certain  given  value.  For  this  reason  the  designer  must 
F-rTS-7327-RE  568 


know  how  to  vary  the  loading  conditions  in  order  to  keep  the 
Maxinun  pressure  constant. 

The  relations  derived  above  permit  one  to  establish  the 
analytical  conditions  under  which  the  maximum  pressure  pB  will 
reaain  constant  when  the  weight  of  the  charge  or  its  density  are 
varied  in  a  given  gun. 

Indeed, 


Pi  fA  F_(9)fu«fq 

—  y  (8)  -  - - ~  - . 

®  1  r  an  .s2I^g 


(61’) 


For  a  given  type  of  powder  (f,  a,  u^)  and  a  given  projectile 
weight  q,  pa  can  be  changed  by  varying  either  A  or  IK  -  e./uj.  A 
is  increased  simultaneously  with  2e^  or  1^,  pR  can  be  kept  constant. 
The  condition  of  the  constancy  of  the  pressure  is' obtained  in  the 


for*: 


P»  - 


»,<•>  w 


const. 


1  -  aA 


Let  us  group  together  the  constants: 


■p^F2(6)  -  const, 


Designating: 
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P — *2^  "  *n  2:5  const' 


we  obtain  the  condition  for  maintaining  p^  constant: 


B  I  —  -  a»a  «  const. 

A  /  m 


(74) 


In  the  case  when  a  /  l/fi,  when  y-  -  0  and  K-  1,  the  condition 

I1 

p  -  const,  has  an  analogous  form:- 


B 


/ 

M- 


const. 


C 


Knowing  and  given  n ,  one  nay  find  the  quantity  B  insuring 
the  obtainment  of  the  given  and,  knowing  B,  one  nay  find  the 
corresponding  value  of  2e^  or  1^,. 

Condition  (-74)  shows  that  in  order  to  keep  pB  constant  when  & 
is  increased,  it  is  necessary  to  increase  the  thickness  2e^  of 
the  powder  in  order  to  offset  the  decrease  of  B  obtained  from 
increasing  the  weight  w  of  the  charge  together  with  the  increase 
of  A  . 

Prom  the  condition  (74)  of  the  constancy  of  the  maximum 
pressure  for. a  giveu  gun,  projectile  and  powder  of  definite  physico¬ 
chemical  properties  (f,  a,  5,  u^),  a  direct  relation  may  be 
established  between  the  weight  **>  of  the  charge,  the  thickness  2e^ 
of  the  powder  or  its  pressure  impulse  1^,  and  the  reduced  length 
of  the  free  space,  in  the  chamber  at  the  end  of  burning: 


■  f0(l  -  QA)  -  *~(1  -  a ^  ) 


1L 

8 
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Indeed,  substituting  the  value  of  B  in  (74)  and  replacing  &  in 
>the  denominator  by  w  /*,  we  obtain: 


s2l|w0(l  -  at)  k2I2W1 

—  mn-T.!!  -  it  r  — ■  i  »  '•  r 


Transposing  all  the  constants  to  the  right  side,  and  designating 


them  by  we.  obtain: 


I-Jw  a.f**  f2F2(6)<em 

— — ' jl  am  .  mm  ■■  "i.  * 


(75) 


Computing  first  from  the  loading  conditions  by  the  following 


■equation: 


- 


F2(6)fV 


and  calculati  g  the  value  «•»  of  the  charge  necessary  to  insure  a  given 
initial  (muzzle)  velocity  vA,  we  can  determine  the  full  pressure 
impulse  1^  «  el/ul 


h  'vV-" 


(76) 


This  equation  shows  that  in  order  that  the  pressure  in  a  given 
gun  remain  constant  when  the  weight  of  the  charge  is  increased,  the 
full  pressure  impulse  Ij,  or  the  thickness  2ei  oi’  powder  must  be 
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increased  at  a  somewhat  higher  rate  than  the  weight  of  t'ne  charge. 

Let  us  apply  equation  (76)  to  calculate  the  powder  thickness 
for  a  76  m w  gur,  1902  model. 

The  conditions  of  loading  are: 


W0  -  1.654;  s  -  0.4693;  0.930;  q  -  6.5;  f  -  9-10',  a  -  1; 

u.  -  7.5-10'6;  1.08;  6  -  0.2;  p  -  2320-102. 

X  £ 


>/r- 

V. 


/o.32*92*1010*l. 08*0. 0663 


*2Pm 


0.4693^*2. 32* 10* 


^36.4*104  -  603; 


ei  n/k^  603-0.930 

K  U1  \/w0  ~  awi  \/l  654-0. 930 


603-0.-930 

0.85 


■  660  kg* sec/dm2; 


2ex  -  2u1Tk  -  2*7.5* 10”6 * 660  -  0.0099  dm  -  0.99  mm  «  1  mm(*) 

We  have  obtained  the  thickness  of  strip  or  tubular  powder  of 
grade  SP  which  was  used  in  this  gun,  and  developed  a  velocity  vA  - 


(*)  This  thickness  of  tubular  powder  corresponds  to  the  grade 
7/7,  while  the  thickness  of  1.28  mm  corresponds  to  grade  9/7. 
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Indeed,  substituting  the  value  of  B  in  (74)  and  replacing  A  in 
the  denominator  by  we  obtain: 


«2lKw0» 

fuifitiii 


2t2„ 

S  *j(*l 


m  i 

r  . 


Transposing  all  the  constants  to  the  right  side,  and  designating 
the*  by  we  obtain: 


1KW1  *****  f2F2(e)<fm 


Computing  first  from  the  loading  conditions  by  the  following 


equation: 


E  - 


?2(fl)f  <fm 


and  calculating  the  value  **»  of  the  charge  necessary  to  insure  a  giveii 
initial  (muzzle)  velocity  vA,  we  can  determine  the  full  pressure 
impulse  -  ci/ui 


1  - 
I 


v 


fl  —  ct  u> 
0 


(76) 


This  equation  shows  that  in  order  that  the  pressure  in  a  given 
gun  remain  constant  when  the  weight  of  the  charge  is  increased,  the 
full  pressure  impulse  1^  or  the  thickness  2e^  of  the  powder  must  be 
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increased  at  a  somewhat  higher  rate  than  the  weight  of  the  charge. 

L  f  us  apply  equation  (76)  to  calculate  the  powder  thickness 
for  a  76  mm  gun,  1902  model. 

The  conditions  of  loading  are: 

Wo  -  1.654;  s  -  0.4693;  0.930;  q  -  6.5;  f  -  9-105,  a  -  1; 

ul  -  7.5*  10~6 ;  1.08;  6-  0.2;  pa  -  2320-102. 


2e^  ■  2u^I^  -  2* 7. 5*  10”^ •  660  «•  0„0099  dm  ■*  0.95  mm  ss  1  nun(*) 

We  have  obtained  the  thickness  of  strip  or  tubular  powder  of 
grade  SP  which  was  used  in  this  gun,  and  developed  a  velocity  vA  - 

(*)  This  thickness  of  tubular  powder  corresponds  to  the  grade. 

7/7,  while  the  thickness  of  1.28  mm  corresponds  to  grade  9/7. 

I 
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>  v>  A  V  v 
'  *  4 

'  ’  ' 

*4  '  ^  n/o- 

t,  '  *) 


»  '588  a/sec  when  the  charge  u>  was  0.930  kg. 

The  same  gun  may  be  fired  with  a  charge  ^  -  1.08  kg  and  a 
velocity  ot^ vQ  620  m/sec  can  be  obtained  at  the  same  pm. 

Find  from  equation  (75)  the  thickness  of  the  powder  in  this 


case. 


Inasmuch  as  \/K^  remains  the  same  as  in  the  first  case, 

■c  *  f 

l 

603*1.08  603*1.08 

r  -  -— =====  -v- - -  -  860, 

VI.  654-1. 08  0.757 

the  thickness  of  tubular  .powder  for  the  same  u^  will  be 

2ex  -  2-7.5* 1G“6* 860  -  0.0129  dm  -  1.29  mm, 

and  this  is  the  thickness  of  our  previous  powder  C^* 

The  ratio 


*K2  _  860  ^  9 

Ijq  660  ’  7  * 


Calculations  show  that- the  equations  derived-  from  approximate 
relations  yield  results  which  are  close  to  experimental  data,  and 
may  be  used  to  calculate  the  variation  in  the  powder  thickness 
concomitant  with  variations  in  the  charge,  if  the  maximum  pressure 
is  kept 'the  same. 
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3.  THE  POSITION  OF  MAXIMUM  PRESSURE  pm  IN  THE  BORE  OF  THE  GUN, 

OR  THE  PATH  /  TRAVERSED  BY  THE  PROJECTILE  AT  THE  INSTANT 

,  HI 

OF  MAXIMUM  PRESSURE, 

,  tu  -  f1  /~F1(e)  -  1 J  -  r0(l  -  aA)/“F1(e>  -  1 J.  (77) 

For  a  given  loading  density  in  a  gun  employing  a  given  type  of 
powder,  the  maximum  pressure  pw  is  developed  at  the  same  distance 

from  the  starting  point  of  the  projectile  (l  -  const.),  regardless 

m 

of  the  powder  thickness  and  the  projectile  weight.  That  is,  the 
path  traversed  by  the  projectile  up  to  the  instant  of  maximum 
pressure  does  not  depend  on  the  powder  thickness  2e^  nor  on 
the  weight  q  of  the  projectile. 

In  a  given  gun  (WQ,  s,  l  Q)  using  a  given  type  of  powder  (f, 

a,  u.,  0)  the  position  (/  )  of  maximum  pressure  depends  only  on 

J.  m 

the  density  of  the  charge.  The  larger  A,  the  smaller  will  be 
(1  -  olA)  and  the  nearer  to  the  starting  point  of  motion  will  be 
p„(*). 

Equation  (77)  shows  that  under  the  condition  of  constancy  of 
the  maximum  pressure  pn,  when  A  is  increased  with  the  simultaneous 
increase  of  B  and  of  the  powder  thickness  2e^,  the  maximum  pn  shifts 
toward  the  origin  of  the  projectile  motion.  Because  the  parameter 
B  increases  thereby,  then,  on  the  basis  of  formula  (62) : 

(*)  This  conclusion  is  also  confirmed  by  the  GAU  tables  for  the 
case  of  Pq  -  300  kg/cm2,  * 1.06,  and  Q  f  i/<5. 
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the  portion  of  the  charge  burned  at  the  instant  of  maximum  pressure 
decreases. 

4.  END  OF  BURNING  AND  PATH  t  TRAVERSED  BY  PROJECTILE 

*  .  K 

< (  IN  THE  BORE  OF  THE  GUN. 

»  The  location  of  the  projectile  at  the  end  of  burning  is  determined 

by  equation  (65),  while  the  corresponding  pressure  p  is  found 
fro*  equation  (66) . 

Fro*  (65)  we  get: 


The  equation  shows  that  for  a  given  loading  density  the 
quantity  increases  with  the  increase  of  parameter  B,  i.e.,  mainly 
with  the  increase  in  powder  thickness.  When  the  weight  of  the 
projectile  is  diminished  while  A  remains  the  sane,  the  path  traversed 
by  the  projectile  at  the  end  of  burning  is  shifted  toward  the 

given  loading  density  the  values 
is  varied  (i.ef,  when  the  powder 
are  varied)  lie  on  the  curve  p,  y 


muzzle  face. 

Equation  (66)  shows  that  for  a 
of  at  the  end  of  burning,  when  B 
thickness  and  the  projectile  weight 
whose  equation  is: 
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This  curye  is  of  the  same  type  as  the  adiabatic  curve  for  . 
instantaneous  burning  whose  exponent  is  however  1  +  6/2.  This 
curve  is  known  as  the  pressure  curve  of  completely  burned  powder. 
The  above  statements  are  clarified  by  the  graph  of  fig.  145. 
The  curves  a,  6  ,  6  and  %  depict  the  pressure  variation  during 
the  burning  of  powders  of  different  thicknesses  when  A  is  the  same. 


Tig.  145  -  Pressure  Curves  for  Different  Powder  Thicknesses. 

a  «  thick  powder;  *>  -  normal  powder;  6  -  thin 
powder;  *-  very  thin  powder. 

2...  A  -  curve  p^y^;  1**.  o  -  curve  p,  y  for 
instantaneous  burning. 

F-TS-7327-RE  -576 


m 


GRAPHIC  NOT  REPRODUCIBLE 


They  are  disposed  in  such  a  way  that  their  maxima  are  on  the 

I 

same  ordinate  at  a  distance  y„  from  the  origin.  The  end  of 

/  *  ' 

burning  occurs  at  a  distance  which  is  governed  by  the  powder 

thickness,  the  'distance  being  the  greater  the  thicker  the  powder 

<y„_  >  rv ,  >  );  the  pressure  values  at  the  end  of  burning 

)Ut’  is 

increase  as  the  powder  thickness  decreases  (p„  <  p_  <  p_.).  The 

points  corresponding  to  the  end  of  burning  lie  on  curve  2-2 

calculated  from  equation  (66)  /"when  6-  0.2,  1  +  0/2  -  1,1  7* 

Curve  1-1  corresponds  to  the  adiabatic  variation  of  the  pressure 

at  instantaneous  powder  burning  (1  +0  «  1.2). 

The  disposition  of  curves  1-1  and  2-2  shows  that  the  pressure 

curves  for  gradual  burning  (a,  $,  a,  %  )  intersect  the  curve  1-1 

depicting  instantaneous  burning.  The  second  period  curves  for  the 

cases  a,  a,  a,  and  t ,*  which  are  not  represented  on  the  diagram, 

are  all  disposed  below  the  curve  2-2  and  above  the  curve  T-l. 

At  the  same  time,  since  for  the  given  powder  p^/(l  -  B8/2)  -  const, 

the  nature  of  the  pressure  change  p  in  the  second  period  depends  upon 

l+O 

the  variation  of  the  variable  factor  1/(1  +  y)  ,  the  latter 
varying  as  in  the  case  of  instantaneous  burning,  i.e.,  along  the 
adiabatic  curve  with  initial  pressure 


- _± -  >  p  . 

M  ' 
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When  the  powder  thickness  is  decreased,  B  and  p^  decrease 
also,  arid  inasmuch  as  the  adiabatic  curves  with  the  same  exponent 
1+8  do  not  cross,  the  adiabatics  in  the  second  period  are 
disposed  the  lower,  the  thinner  the  powder,  i.e.,  inversely 
to  the  disposition  of  the  pressure  curves  in  the  first  period. 

If  we  compare  the  expressions  for  pressure  in  the  second 
period  and  at  instantaneous  burning,  keeping  the  value  of  y  the 
same,  we  will  get  the  following: 

In  the  case  of  gradual  burning  in  the  second  period 


\ 

i 


1 

(1  +  y)1+e 


\ 


ratio  decreases  when  B  decreases  and  the  projectile  weight  increases. 

The  graphs  and  equations  presented  above  for  the  simplest  case 
(#  •»  1,  *Vq  •  0,  a  •  l/6)  permit  one  to  estimate  directly  the 
appearance  and  the  form  of  the  basic  relatipns  between  the  ballistic 
elements  of  a  shot.  They  depict  the  location  and  magnitude  of 
maximum  pressure,  its  dependence  on  the  loading  conditions  (A,  B) 
the  position  of  the  projectile  at  the  end  of  burning  and  the 
gas  pressure  developed  thereby,  the  condition  of  maintaining  the 
maximum  pressure  constant  when  the  weight  of  the  charge  and  the 
powder  thickness  are  varied,  and  the  independence  of  the  curve  y ,/ 
in  the  first  period  of  the  weight  of  the  projectile. 

Such  simple  relations  are  not  obtained  for  the  more  complex 
cases  (*<  /  1,  'Vq  ^  0,  a  /  T/&).  In  such  a  case  it  becomes 

necessary  to  analyze  the  effect  of  the  individual  elements  by 
computing  a  series  of  variations  or  by  using  the  data  found  in 
ballistic  tables. 
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CHAPTER  5  -  A  SURVEY  OF  CERTAIN  OTHER  METHODS  OF  SOLUTION 
(Written  by  Prof.  G.Y.  Oppokov) 

1.  A  VARIATION  OF  PROF.  G.V.  OPPOKOV'S  SOLUTION 
In  order  to  integrate  equation  (13),  Chapter  1,  (p.  473): 


Bxdx 


V  '  +0  +-l!lX  "  B1J 


(78) 


it  is  convenient  to  apply  the  usual  method  of  classical  mathematical 
analysis  -  the  method  of  substitution,  namely,  of  temporarily 
introducing  into  the  process  a  new  variable,  C  *  so  that: 


/  .  c  +  a  &L  x2 

2  ^ 


where  a  is  the  difference  between  the  lengths  of  the  free  volumes 
of  the  chamber  at  the  star tand  end  of  burning: 


When  this  substitution  of  variables  is  effected  in  the  new 
equation,  it  will  become  presently  apparent  that  "the  last  term" 
does  not  contain  in  the  denominator  the  difference: 


%  *  klx  “  Blx  * 


This  obviates  the  need  for  mathematical  transformations  in  the 
course  of  integration  for  the  purpose  of  replacing  the  obtained 
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in'tegral  by  other,  more  simple  ones. 

Indeed,  it  follows  from  equation  (79)  that 


SL  .  +  sea*. 


because ■  /  is  a  constant  in  every  concrete  case.  Moreover: 


+  f  “  ‘(‘Vg  +  kjX  +  *Xx2)  +  r . 


Substituting  in  the  above  the  value  of  /  from'  (79)  : 


L+  7  -  $  +  a  ?|~x2  -  a(4*0  +  kxx  .+  kXx2) 


*  ( 


t  +  l  -  $  -  aty0  +  kxx  -  BjX2) 


We  shall  substitute  into  the  differential  equation  (78)  the 
obtained  values  of  d//dx  and  (l^+  f): 


Bx  • 


+  B0ax  - 


to  +  klx  V 


£\  -  a(t0  +  XjX  -  b^)J. 


We  now  remove  the  brackets  and  effect  the  necessary  simplifications 


+  B0ax  *• 


Vq  +  k^x  -  BjX*\ 


~  K  -  Bax. 
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Grouping  the  verms: 


1L  >  BX 

dX  *Vq  +  kix  “  B  *2 


S  -  -Bad  +  0)x.  * 


*  The  common,  integral  of  this  linear  differential-  equation  of 
th y  first  order  including  the  last  term  can  be  represented  in  the 
following  form: 


'  Bxdx 

WV2 


Bxdx 


?  «  e 


rcx  Ba(l  +e  )  |  e 


’to+kjX-B^B 


xdx  7,  (81) 


where  e  is  the  base  of  natural  logarithms,  and  is  still  an  arbitrary 


constant. 


We  must  introduce  into  the  analysis  Prof.  Drozdov’s  function: 
C  Bxdx 


'fyj+kjX-B^X'' 


Z 


Then  the  partial  integral  of  the  last  equation  with  respect  to 
the  derivative  dS/dx  will  be: 


B1  f  B1 

S  -  Z  £'Cl  -  Ba(l  +6)  \  Z 


xdx  /. 
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We  shall  now  return  to  the  desired  path  l  of  the  projectile, 


t 

\ 


C 


( 


and  after  substituting  the  obtained  value  of  £ ,,  get: 


(  -  Z 


B 


B 

s: 


Ba  (1  +  6 ) 


,  —♦  B  8  2  » 

xdx  /  +  a  —  x  -  f 
2 


Let  us  determine  now  the  constant  C,  from  the  initial  conditions 

1 


at  which: 


B 


t  -  0;  x  -  0;  Z 


1  f 

-  i; 


Zxdx  «  0, 


We  have  from  the  latter  equation  for  the  path  of  the  projectile: 


1  G> 


0  -  1  (C  +0)+0-|, 

1  A 


whence 


C1  "  lb' 


Thus  the  desired  path  of  the  projectile  is  defined  by  the 
following  expression: 


JL  * 
B1 

z  ft 


x 


B 

s: 


-  aB(l  +  6) 


,  —f  B0  2  f 

xdx  /+  a  — — x  - /  . 
—  2  ^ 


(82) 
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2.  PARTICULARS  OF  PROF.  I.P.  GRAVE  'S  SOLUTION' 


This  method  of  solution  was  developed  by  Prof.  I.P.  Grave  in 
order  to  perfect  Bianchi's  method,  which  was  the  first  variant 
(in  time)  of  the  / ^  method.  In  Bianchi's  original  equations  the 


av. 


effect  of  the  variation  of  l 


is  discounted  and  in  integrating 


av. 


he  considers  the  quantity  as  a  certain  incompletely  determined 

constant.  Bianchi  divides  the  curve  p,  l  into  three  segments, 

for  which  instead  of  he  takes  ~a  and  l  y  respectively, 

*  2 

which  corresponds  to  the  following  conditions: 


H*  -  0;  'V  «  0.5;  «  1. 

av.  av.  av. 


But  in  this  case  the  curves  p,  I  and  y,  [  are  not  smooth;  they 
have  angular  points  corresponding  to  the  beginning  of  the  second  an a 
third  segments. 

In  order  to  take  into  account  the  effect  of  the  variation  of 
and  obtain  smooth  p,  {  and  v,  /'  curves,  Prof.  I.P.  Grave,  in 
integrating  the  equation  for  the  path  of  the  projectile,  considers 
as  a  variable,  defined  by  the  average  value  of  its  derivative  with 
respect  to  l  .  This  average  value  of  the  derivative  must  be  negative, 
because  decreases  during  the  burning  of  the  powder.  Consequently: 


-k, 


whence,  after  integration,  we  obtain: 


i 


Hr 


-k  ( . 


(83) 
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and,  consequently,  the  following  must  obtains 


;  -(l-k)e- 

(i  -  k)  i-  - 1*  z  - 1  . 

Using  a  second  time  equation  (83),  we  obtain  the  following 
equation  from  the  above: 

-(l-k)g- 

l"  L  Z  -  /*•  (85) 

A  certain  difficulty  arises  from  the  fact  that  in  order  to 

apply  equation  (85)  it  is  necessary  to  know  the  constant  k,  for 

which,  in  turn,  it  is  necessary  to  know  1%/lj,  /"see  equation  (84)_7; 

but  the  path  („  of  the  projectile  at  the  end  of  the  period  is 

It  •- 

unknown  beforehand. 

In  order  to  overcome  this  difficulty,  a  nomograph  is  given  • 

/~ see  I.P.  Grave,  "VNUTRENNYAYA  BALL I ST I KA"  (Internal  Ballistics), 
Pyrodynamics,  No.  1,  p.  58_7  which  enables  one  to  determine  the 
ratio  I if  * - 

l+n  $(xK)  B  ~1  B  l 

and  - —  -  -—In  ZY  -  2.303  —  log  Z~  . 

<1  H#X  Bl  *  K 

are  known. 

Having  found  from  the  conditions  of  loading:  * 

,  •  ! 
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~t —  and  2 .  J03  « —  lo^j  Z  , 

ll  1  K 


we  can  determine  1  j^/  /  ^  from  this  graph;  this  will  enable  us  to 
calculate  k  from  equation  (84). 

If  it  is  necessary  to  find  k  raore  accurately,  the  obtained 
value  of  k  may  be  rendered  more  exact  by  successive  approximations. 
The  value  k^  found  from  the  graph  is  substituted  into  (85), 
is  found  in  the  second  approximation,  and  a  new  value  k2  is  then 
determined  from  (84)  representing  a  second  approximation,  and  so  on, 
until  two  consecutive  values  of  coinciding  with  the  required 

degree  of  accuracy  are  obtained. 
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CHAPTER  VI  -  SOLUTION  OF  THE  FUNDAMENTAL  PROBLEM  OF  INTERNAL 
BALLISTICS  ON  THE  BASIS  OP  THE  PHYSICAL  LAW  OF  BURNING 

(M.Ye.  Serebriakov’s  Method) 

I 

As  was  shown  in  Part  I  of  this  text,  the*  actual  burning  of 
powders  deviates  from  the  geometric  law  under  the  influence  of  a 
number  of  factors.  An  analysis  obtained  by  the  aid  of  the  progressivity 
curves  P,  ^  and  P,  t  has  shown  that  certain  anomalies  actually  occur 
^even  during  the  burning  of  powders  of  simple  shapes:  non-instantaneous 
ignition,' accelerated  burning  of  the  outside  layers  (ballooning), 
etc.  The  burning  law  cannot  be  established  at  all  on  the  basis  of 
the  geometric  law  for  adulterated  and  porous  powders  used  in  pistol 
cartridges. 

The  actual  burning  law  can  be  established  only  by  burning  powder 
in  a  test  bomb  at  different  loading  densities  and  by  obtaining  pressure¬ 
time  curves  reflecting  all  the  deviations  and  peculiarities  of  a  given 
sample. 

t 

-  The  variation  in  the  intensity  of  gas  formation  T  and  1  -  ^  pdt 

0 

as  a  function  of  vf  and  t  can  be  established  from  the  obtained  p,  t 


curve. 


Both  graphs  p>  H*  and 


^  pdt,  in 


conjunction  with  the  fundamental 


p,  t  curye  obtained  from  the  bomb  test  enable  us  to  soive  the 
fundamental  problem  of  pyrodynamics,  i.e.,  to  compute  the  gas-pressure 
and  projectile  velocity  variation  curves  under  conditions  of  actual 
burning  of  powder  in  the  bore  of  a  gun. 


These  graphs  also  enable  us  to  establish  the  individual  behavior 


of  powder  lots  of  different  grades  occasionally  differing  considerably 
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as  to  their  properties,  which  behavior  could  not  be  disclosed  by  any 
method  other  than  by  bomb  tests. 

There  have  been  actual  cases  where  powder  lots  of  the  same 

grade  and  the  same  manufacturer  having  identical  chemical  composition 

and  dimensions  produced  a  difference  of  6-8%  in  the  charging  weights 

when  fired  at  the  same  values  of  v .  and  p 

^  max 

Bomb  tests  had  shown  that  the  burning  rate  u^  of  these  powders 
varied  as  much  as  15-20%;  This  variation  could  not  have  been  disclosed 
by  any  other  means  except  the  .bomb  test. 

The  fundamental  problem  cf  internal  ballistics  for  adulterated 
or  porous  powders  can  be  solved  in  exactly  the  same  manner  only  on 
the  basis  of  the  experimental  (physical)  law  of  burning. 

We  are  presenting  below  the  method  of  solving  the  fundamental 
problem  on  the  basis  of  the  physical  law  of  burning,  when  the  burning 
rate  law  is  u  -  u^p,  which  corresponds  to  the  coincidence  of  the  curves 
I,  t  or  j  pdt,  ^  at  various  loading  densities. 

The  basic  assumption  made  here  is  that  both  in  a  bomb  at  different 
loading  densities  A  and  in  a  weapon  with  a  variable  space  in  the  case 
of  a  continuously  decreasing  loading  density,  the  value  of  J  pdt 
■  is  a  single-valued  function  of  only,  and  does  not  depend  on  the 
loading  density.  This  condition,  which  has  been  proved  by  bomb 
tests  at  different  loading  densities,  is  being  extrapolated  in  the 
given  case  for  considerably  higher  values  of  A  in  a  weapon. 

1.  DERIVATION  OF  BASIC  RELATIONSHIPS  AS  APPLIED  TO  THE 
PHYSICAL  LAW  OF  BURNING. 

The  solution  is  based  on  applying  the  pressure  curve  obtained 
from  bomb  tests  to  the  computation  of  curves  depicting  the  gas 
pressure  and  velocity  of  the  projectile  in  a  weapon. 
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As  the  projectile  moves  through  the  bore  and  the  initial 
air  space  becomes  larger,  the  pressure  will  depend  on  the  current 
value  of  the  loading  density  A  -•*  +  s tj. 

We  shall  introduce  the  following  designations: 

A  ^  ~  loading  density  of  powder  in  test  bomb,; 

AQ  -  initial  loading  density  in  weapon. 

P  and  T  -  gas  pressure  and  time  corresponding  to  the  given  value 
of  at  constant  loading  density  A^,  at  which  the  bomb 
test  was  conducted  and  at  which  the  curve  P, t  was 
obtained ; 

p  and  t  -  gas  pressure  and  time  corresponding  to  the  same  value 
of  when  A  is  variable,  which  condition  applies  to 
a  given  disposition  of  the  projectile  in  the  bore  of 
the  barrel. 

We  shall  designate  the  corresponding  integral  values  as  follows: 


a)  In  bomb 
* 


b)  In  weapon 


Pdt  -  I 


pdt  ■*  i 


Pdt  -  I, 


pdt  -  i, 


Pdt  -  I, 


pdt  -  iK. 


I  is  obtained  from  a  table  or  graph  as  a  function  of  y  or  t , 
on  the  basis  of  bomb  tests. 
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-  *gr 


We  shall  consider  henceforth  the  value  of  as  the  independent 


variable.- 


2.  DETERMINING  THE  PROJECTILE  VELOCITY  AS  THE  FUNCTION  OF  vy 
•  * 

On  the  basis  of  the  impulse  theorem 
tyradv  -  spdt. 

Integrating  from  the  start  of  motion: 


<pmv  -  s  f  pdt  -  s(i  -  iQ)., 


where  vyQ  is  the  portion  of  the  charge  burned  in  the  gun  at  the  start 


of  notion: 


2_  _L 


—  +  a  -  -~ 


Determining  v: 


v  -  -2-(i  -  i  ) 
<pm  0 


or  on  the  basis  (86) 


is  | 


—  a  -  i) 

«fm  0 


The  values  of  I  and  I  for  vP  and  vy  are  known  from  the  bomb 

0  1  0 

test,  s/«fn  is  known  from  the  gun  data.  Thus  the  velocity  of  the 
projectile  is  determined  from  equation  (89)  as  the  function  ofvy. 

At  the  end  of  burning  when  ^  -  1 
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(90) 


v 

K 


In  contrast  to  the  analogous  formula  in  the  case  of  the  geometric 

law  of  burning,  the  value  of  I  corresponds  not  to  the  average 

K  / 

thickness  of  the  powder  but  to  the  maximum  thickness,  which  may 

considerably  exceed  the  average  thickness  of  powders  having  a  variable 

thickness.  A  diagram  of  the  pressure  impulse  of  tabular  powder 

usually  obtained  in  bomb  tests  ifj  offered  in  fig.  146. 

The  value  of  the  impulse  1,  »  e,  /u  corresponds  to  the 

i  X  AVo  1 

burning  of  powder  of  average  thickness;  the  value  of  I  I 

A 

corresponds  to  the  burning  of  the  thickest  element  of  the  charge. 


r: 


Fig.  146  —  Pressure  Impulse  of  Tubular  Powder 

\ 

Therefore*  also  v  in  the  case  of  actual  powder  burning  assumes 

K  \ 

a  considerably  greater  value  than  in  the  case  of  the  geometric  law, 

.  \ 

t 

but  in  that  case  the  projectile  will  also  traverse  a  considerably 
longer  path  at  the  end  of  powder  burning,  so  that: 


I 


K 


K 


(Physical  law  of  burning) 
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3.  DETERMINING  THE  PATH  OF  THE  PROJECTILE  AS  A  FUNCTION  OF 
In  the  wetbod  outlined  the  relation  between  t  and  <y  is  established 
by  means  of  the  auxiliary  function  L,  14/  determined  from  the  same  bomb 
test  at  a  loading  density  A^i  by  the  additional  analysis  of  the 
test  curve  P,  T  .  And  if  the  value  of  the  pressure  Impulse  I,vy 
does  not  depend  on  the  choice  of  A^,  then  the  function  L,  vy  depends 
on  the  value  of  Ai  chosen  at  the  test. 

As  we  shall  see  later,  the  function  L  has  the  dimensionality 
of  the  path  and  a  definite  physical  meaning. 

We  will  have  from  equation  (88)  : 


d{  -  vdt  -  --£-(i  -  ift)dt,. 


where  (i  -  i0)  is  a  function  ofvy. 

Upon  integrating,  the  path  of  the  projectile  will  be  determined 
by  the  formula 


* 


(91) 


Here  the  element  dt  corresponds  to  the  element  d*  when  the  powder 
is  burned  under  conditions  of  variable  volume  (space)  and  depends  on 
the  value  of  pressure  p  at  any  given  instant  which  is  still  unknown. 

We  can  obtain  an  expression  from  the  bomb  test  analogous  to 
expression  (91) : 
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s 

m 


G. 


* 


02) 


The  value  of  L  *•  function  of  y  -  is  obtained  by  the  second 

¥ 

integration  of  curve  ^  Pdt  with  respect  to  *c  and  by  multiplying 
same  by  the  coefficien?  s/^m.. 

If  t  is  the  path  traversed  by  the  projectile  at  the  instant 
the  portion  vy  of  the  charge  is  burned,  then  L  is  the  path  the 
projectile  would  have  traversed  if  the  pressure  behind  it  developed 
according  to  the  same  law  as  in  a  boob  with  a  constant  loading 
density  at  the  instant  the  same  portion  of  the  charge  vy  is 
burned. 

«* 

L  has  a  definite  physical  meaning.  For  example,  it  is  obtained 

in  practice  in  a  bomb  of  considerable  capacity  with  a  free  piston. 

3 

Thus  in  a  bomb  of  capacity  -  300  cm  the  piston  displacement  (with 

2 

the  piston  usually  having  a  cross-sectional  area  of  s  -  1  cm  )  is 
about  3  c*.  The  change  in  volume  amounts  to  only  1%,  and  hence  the 
piston  is  displaced  by  a  pressure  which  increases  in  almost  a  constant 
volume. 

The  value  of  L  as  a  function  is  found  from  the  bomb  test  using 
the  procedure  given  in  the  table  below. 

It  is  necessary  to  establish  the  relation  between  {  and  L, 
and  hence  between  (  and  because  L  is  function  of  ^ . 

Differentiating  equations  (91)  and  (92)  and  taking  their  ratio 
and  reducing,  we  get: 
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(63) 


From  equation  (87) 


whereas  the  ratio  P/p  is  replaced  by  the  ratio  of  the  free  volumes 
on  the  basis  of  the  equation  of  state: 


RT,u> 


1*1* 


fAj^ 


—  -  A,  a 

&  1  V 


lb 


T  -  “ 


The  expression  in  the  denominator  represents  the  free  specific 
gas  volume  at  loading  density  A^. 

An  analogous  expression  will  obtain  also  in  the  formula  for  p. 
We  shall  replace  the.m  with  average  values,  because  the  last  term 
is  small  in  comparison  with  th >  first  two,  and,  moreover,  the  ratios 
of  the  free  volumes  will  enter  hem  all. 


-  Wn  1 - 1  _  a  (  a  -  -L  L 

av.  0  8  1\  5  /Tav. 


-  M 1  ~ 


- A 

2  1 


-  W0(l  -  a’Ax), 


where  ex*  » 


1 

a  +  T 
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For  a  gradually  burned  powder 


*o  -  tt,“ 

if  -  a '  w  +  si 
0 


Vo 


The  ratio  P/p,  following  substitution  of  the  proper  expressions 
and  simplification,  will  take  on  the  form: 


Upon  incorporating  this  expression  in  (94)  and  then  in  (93), 


we  ge  t : 


*0  (  t  \  © 

dl  -  dL  -  {  1  +  — —  I  ,  where  k*  «  1  +  — . 


Dividing  the  variables  and  integrating: 


-K  r  'M 
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Designating  1  -f  W  “  x,  we  get  in  the  left  sido^r 


l 


i 

arid  in  the  right  side 


We  then  obtain: 


0  but 

r 


U 
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6  S 
~2  “ 


-  a’ 


Solving  (96)  for  J,  we  get: 


JL 

B'L'*’  )  e 
■*0 


As  in  the  solution  of  the  problem  of  internal  ballistics  by 
the  average  |  method,  the  value  l  -  average  for  the  entire  burning 

iT  C 

‘process  -  can  be  replaced  in  this  formula  by  the  current  value  / 

't’av. 

by  means  of  the  usual  formula: 


tQ  l  _  A.  -  A  (a  -  -- 

0  ft  \  A 


vy  4.  vy 


and  then 
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‘  *  At  jJt 


Formula  C100)  gives  the  path  I  as  a  function  of  y  by  means 
of  the  auxiliary  function  L*  ,  determined  in  bomb  tests  at  loading 

density  A  «  which  function  reflects  (depicts)  ^he  true  burning 

law. 

/ 

Co. paring  this  formula  with  the  analogous  formula  used  in  the 
method  to  solution  in  which  I  ,  we  will  „oto  that  in  place 

of  Prof.  Drozdovs  function  Z%  1  formula  (100)  contains  the  function 


(1  -  B'L*  ) 


-2/e 


To 


For  the  case  where  X  -  1,  X-  0, 


B 


Bt 


_2_ 

e 


Therefore  the  expression  in  parentheses  (1  -  B'Lj  )  has  replaced 

in  this  solution  the  function  Z  in  the  case  of  the  geometric  law 
•  * 

{  of  burning. 

v 

Formulas  (89)  and  (100)  enable  us  to  calculate  and  plot  the 
projectile  velocity  curve  as  a  function  of  path  l  . 

Pressure  p  is  found  from  the  fundamental  equation  of  pyrodynamics : 

(101) 


f  e  2 
fwt - g  Cfmv 


P  “ 


8(^  +  t  ) 


wherein  the  variable  quantities  as  the  y  functions  are  already  known. 
To  determine  the  maximum  pressure  p*  and  the  corresponding  value 
C  *e  differentiate  the  equation  (101)  with  respect  to  t: 
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the  geometric  law,  we  will  obtain  the  usual  relationship  for  x  . 

» 

At  the  end  of  burning  when  ^  -  1,  we  will  have: 


A 
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Ski- 


Fig..  147  -  Determining  for  Maximum  Pressure. 


4.  GRAPHICAL  CLARIFICATION  OF  THE  METHOD  OF  SOLUTION 

In.  order  to  solve  the  problem  on  the  basis  of  the  physical  law 

of  burning,  it  is  first  necessary  to  perform  the  ballistic  analysis 

of  the  given  powder.  To  do  so,  bomb  tests  are  conducted  at  two 

loading  densities  A  and  A  >  the  ballistic  characteristics  - 

1  2 

propellent  force  of  powder  f  and  covolume  a-  -  are  determined,  and 

also  the  test  characteristic  of  the  intensity  of  gas  formation 

t 

and  the  impulse  of  pressure  increase  f  PdT  -  I  (fig.  148). 


Knowing  the  loading  conditions  of  the  weapon,  we  determine 
1  1 


*0  “ 


■;  from  graph  I,  we  find  the  corresponding  value 


t  1 

-  +  a - - 

p0  i 


Otlo-  j 


PdT,  and  upon  subtracting  this  value  of  1^  from  all  the 


values  of  I,  obtain  the  dependence,  of  I  -  IQ  on  ^  and  T  (fig.  149). 
Integrating  numerically  the  curve  I  -  I  ,  T  with  respect  to  X ,  we 

y  •  0 

find  the  integral  ^  (I  -  IQ)df  as  a  ,  function  of  . 

*0 
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Introducing  the  designation: 


i  f  PdtdT  -f(l-  v 


Multiplying  I  -  In  and  G*  by  s/<fm,  we  get 

u  % 


V  -  _  (I  -  1  ); 

re  m  V 


(89) 


-  G'f  - 

<pm  *0 


8  * 

*0 


a  -  Vdt 


(92) 


▼  does  not  depend  on  A (a^  or  a 2>  and  is  a  function  of  y  only 
for  all  the  loading  densities.  The  function  ,  being  a  function  of 

depends  at  the  same  time  on  A,  because  the  time  element  dl  during 
which  a  definite  portion  of  charge  dvy  is  burned  decreases  with  the 
increase  of  A . 

Indeed,  from  the  equality 


it  follows  that 
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j_.  dY 

dr- 

rp 


A.  J— i  -A. 


where 


/  i  \  1  - 


and  hence  dt  varies  inversely  with  the  change  of  A,.  / 

1  , 

For  this  reason  the  curves  G1*  and  L*  as  a  function  of  wi'll 
also  be  disposed  the  lower  the  greater  when  testing  the  powder 
in  a  bomb  (see  fig.  149). 


Fig.  148  -  Basic  Curves 
and  I,  vp  . 


Fig.  149  -  Auxiliary  Curves 
for  Determining  the  Elements 
of  a  Shot. 


The . arrangement  of  a  table  for  analyzing  bomb  tests  as  a  means 
of  obtaining  all  the  auxiliary  functions  is  presented  below;  this 
table  serves  to  clarify  the  graphs  in  figs.  148  and  149. 
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Procedure  for  Analyzing  a  Bomb  Test  at  Loading  Density  ^  as  a  Keans  of 

Obtaining  the  Basic  Functions 


We  shall  construct,  recording  to  equations  (89)  and  (92)  for 

the  same  values  of  vy,  a  curve  showing  the  dependence  of  v  on  L  when 

A «  /~fig.  150,  curve  v,  L(A^)_7-  If  the  test  were  analyzed 

for  A  >  A, ,  the  relationship  v,  L(A  )  would  obtain,  which  relationship 

curve  has  a  larger  slope  angle  and  in  which  the  shorter  path  of  the 

projectile:  L„(A  )  <  L  (A  )  corresponds  to  the  end  of  powder 
xv  2  K  i 

burning.  Both  curves  have  a  smaller  curvature  than  the  true  v,| 

curve  of  the  velocities  of  the  projectile  in  the  bore. 

If  we  were  to  extrapolate  the  function  L^,  for  the  initial 

0 

loading  density  a  in  the  gun,  we  would  have  obtained  curve  v,  L(a  ) 

0  0 
(150),  which  at  the  start  of  motion  has  a  common  point  of  tangency 

with  the  true  v, l  curve.  The  latter  is  obtained  when  A  decreases 

continuously,  and  as  and  v  increase  curve  v,  l  (heavy  dotted  line) 

gradually  goes  over  from  curve  v,  L(AQ)  to  curves  corresponding  to 

ever  smaller  a,  which  family  of  curves  includes  also  curves  v,  L(a2) 

and  v,  L(ax) . 

This  transition  is  the  one  given  by  the  fundamental  formula 
(100): 

2 

i- 1  ra  -  b'ia  >"  T  - 1  7 

Tav.~  *0 


together  with  formula  (89) : 


cpm 


d  -  y 


The  difference  between  the  method  of  solving  the  fundamental 
problem  outlined  above  and  other  such  methods  lies  in  the  fact  that 

t 
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ia  deriving  the  dependence  of  l  on  \y,  use  is  made  not  of  the 
fundamental  equation  of  pyrodynaraics,  but,  rather,  of  the  equation 
of  the  state  of  powder  gases  for  different  positions  of  the 
projectile  in  the  bore  of  the  barrel. 

In  solving  this  problem  use  is  made  of  the  gas  pressure  curve 
P,  X  obtained  from  bomb  tests,  which  expresses  the  true  burning  law 
with  all  the  deviations  from  the  geometric  law. 

An  analogous  result  can  be  obtained  only  by  the  numerical 
integration  of  Taylor's  series  or  from  finite  differences. 


Fig.  150  -  Relation  Between  Auxiliary  Curves  v,  L  and 

Actual  Curve  v,  /  . 

The  method  outlined  here  permits  the  solution  of  the  problem 
also  in  the  case  of  the  geometric  law,  by  assuming  the  .following 
theoretical  relationship  forr,y: 


Therefore,  this  method  is  a  more  general  one  than  the  methods 
based  on  the  geometric  law  of  burning. (♦) 

(*)  For  a  more  detailed  explanation  see:  Serebriakov,  K. Ye.  "FIZICHESKY 
ZAXON  CORENIA  YO  VNUTRENNEY  BALLISTIKE"  (The  Physical  Law  of  Burning  in 
Internal  Ballistics).  "030R0NGIZ"  (State  Publishers  of  Defense  Literature) 
1940. 
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5.  ANALYSIS  OF  THE  OBTAINED  CURVES  p ,  l  AND  v,  l  . 
Analysis  of  curves  p,  i  ^nd  v,  l  obtained  on  the  basis  of  the 
physical  law  of  burning  indicates  that: 


Fig.  151  -  Curves  p,  /  and  v, /  Obtained  on  the  Basis  of  the 
Physical  and  Geometric  Laws  of  Burning. 

$.5'%.  -  physical  law  of  burning. 

.  3  .  1 .  -  geometric  law  of  burning. 

1)  Due  to  ballooning  -  accelerated  burning  of  outer  layers  - 
the  maximum  pressure  is  attained  earlier,  and  the  pressure  curve  is 
disposed  higher  at  the  start  than  in  the  case  of  the  geometric  law; 

2)  The  beginning  and  the  first  half  of  the  velocity  curve  v,  / 
obtained  on  the  basis  of  the  physical  law,  are  disposed  above  the 
corresponding  v,  {  curve  in  the  case  of  the  geometric  law  of  burning; 
the  curves  merge  at  the  end; 

3)  The  same  value  of  p  is  obtained  at  a  smaller  propellant 

max . 

force  of  powder  than  in  the  case  of  the  geometric  law; 

4)  Due  to  after-burning  of  the  thicker  elements  of  the  charge, 
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the  end  of  burning  is  transposed  nearer  the  muzzle  face,  and  the 

velocity  exceeds  the  theoretical  value,  for  the  average  powder 
thickness;  / 

5)  Due  t6  the  gradual  decrease  of  the  intensity  of  gas  formation 
F,  Mr*  at  the  end  of  burning,  the  transition  of  the  pressure  curve 
from  the  first  period  to  the  second  proceeds  without  d  jump  and 
forms  no  turning  points  on  the  curve,  as  it  does  in  t^ie  case  of  the 
geometric  law,  ' 

The  graph  in  fig.  151  clarifies  the  above. 
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